VBEHLJPS - 25/N0V/93 



GENERALIZED HARTREE-FOCK THEORY 
AND THE HUBBARD MODEL 



Volker Bach^^), Elliott H. Lieb(2)'(3) and Jan Philip Solovej^^) 

_ W FB Mathematik, Technische Universitdt Berlin, Strasse des 17 Juni 136, D- 10623 Berlin, Ger- 

m . ' ' ' ' 

Q"^ , many 

0^ \ ^"^^ Department of Physics, Jadwin Hall, Princeton University, P.O. Box 708, Princeton, NJ, 08544 

^^"^ Department of Mathematics, Fine Hall, Princeton University, Princeton, NJ, 08544 

Q 

o ■ 

Dedicated to Philippe Choquard on his sixty fifth birthday 

> 

, Abstract: The familiar unrestricted Hartree-Fock variational principle is generalized 

to include quasi-free states. As we show, these are in one-to-one correspondence with 
the one-particle density matrices and these, in turn provide a convenient formulation 

■ of a generalized Hartree-Fock variational principle, which includes the BCS theory as a 
special case. While this generalization is not new, it is not well known and we begin 

! by elucidating it. The Hubbard model, with its particle- hole symmetry, is well suited 

to exploring this theory because BCS states for the attractive model turn into usual HF 

■ states for the repulsive model. We rigorously determine the true, unrestricted minimizers 
C , for zero and for nonzero temperature in several cases, notably the half-filled band. For 
Q I the cases treated here, we can exactly determine all broken and unbroken spatial and 

gauge symmetries of the Hamiltonian. 
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I. INTRODUCTION 

The Hubbard model has long been recognized as an interesting imitation of electron-electron 

interactions and of the correlations they induce. It also turns out, as we show here, that it is an 
interesting testing ground for Hartree-Fock (HF) theory, and it is one of the very few examples in 
many-body theory for which many properties of the true, energy minimizing HF ground state and 
pressure maximizing positive temperature state can be precisely elucidated without approximation, 
restriction or unjustified assumptions. (To avoid misunderstandings, we should make it clear at 
the outset that we consider only unrestricted HF theory.) 

While studying the HF theory of the Hubbard model we were led to a critical study of HF 
theory itself — namely, the proper context in which to view it, as well as some of its very general 
features. This is the content of Sect. II which we can summarize as follows. 

The usual HF theory for an A'^-particle system starts with a Slater determinant -0 formed 
from N orthonormal, one-particle functions of space and spin; the energy is then minimized with 
respect to the choice of these N functions. This ip will generally break certain symmetries inherent 
in the original problem — typical examples being translation invariance and spin or spatial angular 
momenta. That being the case, it is not worse to permit a violation of particle-number conservation 
as well — if the energy can thereby be lowered. This is precisely what was done by Bardeen, 
Cooper and Schrieffer [BCS] in their theory of superconductivity. It is important, here to recall 
the simple fact that if a state violates a symmetry of the Hamiltonian (e.g. angular momentum. 
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particle number) and if that state is then decomposed into its irreducible components, at least one 
component will have an energy no greater than (and often less than) the original state. 

The proper Hilbert space, then, is Fock space, the sum of all the original A^-particle spaces 

with iV = 0, 1, 2, The simplest Hamiltonians on T are, of course, the quadratic ones, and they 

are all diagonalizable by Bogoliubov unitary transformations, which transforms creation operators 
into linear combinations of c's and c^'s. In the context of fermionic theories, this transformation 
was discovered at the same time by Valatin [VJ] and Bogoliubov (see [BN]). The ground state of a 
quadratic Hamiltonian H'^ is a Bogoliubov transformation lA applied to the vacuum (zero-particle 
state) |0). We call all such states of the form U\^) a generalized HF state because the original 
A''-particle determinant (which has the form cjcl . . . c]y|0) and which we call a normal state) is 
only the special case corresponding to an H'^ that contains terms of the form c^c and terms 
of the form cc^, but not the particle nonconserving terms cc or c^c^ . Such generalized HF states, 
U\0), are also called quasi- free states because they satisfy the conclusion of Wick's Theorem. 
Indeed, all quasi-free states are of this form, as we show below. 

A generalized HF ground state is thus the ground state of some quadratic in Fock space. 
A positive temperature HF state is, likewise, the usual grand canonical Gibbs state for such an 
H^. In analogy with the ground state, such a grand canonical Gibbs state is called normal when 
H'-^ is particle conserving. For the ground state, the H'^ is determined so that the expectation 
value of the original many-body Hamiltonian, H, of our system is as small as possible. For positive 
temperatures, is chosen (in a temperature dependent way) to maximize the pressure in the 
grand canonical ensemble. 

The paired (BCS) state, which is so important in superconductivity theory, is also of the form 
U\0), a Bogoliubov transformation of the vacuum. Thus, normal HF theory and BCS theory are 
but two aspects of the same general theory: find the best substitute quadratic Hamiltonian or, 
equivalently, find the best Bogoliubov transformation. This relationship was certainly known [BN], 
[dG], [BR], [DJK], but our personal experience is that it is far from being universally appreciated. 
At first it seems surprising that one-particle states can somehow evolve into pair states, and the 
explanation is roughly the following: It is always true that for each mode a, Uc]jA^ = + 
where d) (resp. e) is proportional to a creation (resp. annihilation) operator, but it is possible 
that = or e = 0. If = then ijj = IA\^) contains the mode e^, i.e., e)ei]j = (pos. const)^. If 
e = 0, ip just contains a factor proportional to |0), i.e., dij^ = 0. If e and are both nonvanishing 
then such modes must come in pairs (as we prove in Theorem 2.2) and ip is found to contain a 
pairing factor (l + d^e^) acting on |0). Alternatively, it turns out that to every eigenvalue 1 of the 
particle conserving part, 7, of the one-particle density matrix (1-pdm), T, associated with the state 
ip =U\0) there corresponds a simple one-particle state in tp, but to every eigenvalue of 7 between 
and 1 there corresponds a pair state in ip. 

The 1-pdm plays an essential structural mathematical role in HF theory. The set of quasi- 
free states (generalized or not) does not have a linear or convex structure. A linear or convex 
combination of such states is not necessarily a qTiasi-frce state. However, a convex combination 
of 1-pdms is a 1-pdm, i.e., F = AFi + (1 — A)r2 is a 1-pdm if < A < 1 and Fi and r2 are 
1-pdms. This fact allows us to compensate for the missing convex structure of quasi-free states. 
Thus given -01^2 = ZYi,2|0) with Ui^2 being two Bogoliubov transformations, we can form the 1-pdms 
Fi and r2 and then form F as above. Finally, we can return to the level of the quasi-free states 
and thereby define a quasi-free state that interpolates between the two original states. Sect. II 
contains a detailed description of quasi-free states, density matrices and quadratic Hamiltonians. 
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We present this partly for the reader's convenience, but also because we could not find quite what 
we need in the literature on quasi-free states (which usually concentrates on quasi-free states in 
terms of algebraic automorphisms rather than operators) or in the excellent book by Blaizot and 
Ripka [BR] which does not deal explicitly with the infinite dimensional Hilbert space L^(R'^), the 
space of square integrable functions on R^, needed for our other theorems in Sect. II on atomic 
HF theory. 

The notation in Sect. II is a bit complicated and one reason for this is the necessity to introduce 
antiunitary transformations (because c's transform by antiunitaries if c^'s transform by unitaries). 
Consequently, if one tries to write equations in a basis independent way, one needs more than the 
usual notation of linear operator theory. If one fixes a basis, however, one can use the ordinary 
linear operator notation, but then complex conjugates (denoted by superscript bars) and indices 
appear in profusion. We have opted for the second route. 

The general Theorem 2.14 in Sect. II about the usual (A'^-particle) HF theory should have been 
well known but seemingly was not. It applies to repulsive two-body potentials (as in the real world 
of electrons with Coulomb interaction) and states two things. The first is that the N one-particle 
states are precisely the energetically lowest eigenvectors of the HF operator. (This fact was stated 
in [LS], and the proof was sketched in [LE4]). While the N HF orbitals are distinct eigenvectors 
of the HF operator, it is not obvious, a priori, that they are the lowest ones; indeed, this might 
not be true when the interactions are attractive. The second part is surprising, for its conclusion 
runs counter to what one might naively assume. There are never unfilled shells (for any choice of 
A'^). That is to say, the degeneracy of the last level of the HF operator is always precisely what is 
needed to accommodate the available number of electrons — not more than that! 

Readers who are already familiar with the formalities of Sect. II are advised to skip over it and 
to turn to Sect. Ill, which contains the HF analysis of the attractive Hubbard model. It beautifully 
illustrates the relationship between BCS and normal states discussed above. It is well known that 
a particle-hole transformation (on the down spins alone) converts the repulsive and the attractive 
models into each other. What does this do to generalized HF states? The answer, simply, is that a 
normal HF state may be turned into a BCS state (in which there is pure pairing without isolated 
one-particle states). Indeed, it turns out that the repulsive Hubbard model at half-filling (i.e., the 
expected particle number equals |A|, the number of sites in A) always has a normal state as its 
optimal state (for zero and for positive temperature). The attractive model at half-filling then 
has a BCS state (and, when the lattice is bipartite, also a normal state of the same energy) as 
its ground state. This was well understood by Dichtel et. al. [DJK], but they did not prove that 
their state was, indeed, the true minimum energy HF state. We do so here as a special case of our 
results in Sects. Ill and IV. 

At the outset we emphasize that translation invariance is not assumed. By the word "lattice" 
we mean a collection of points connected by bonds (or edges). Perhaps "graph" would be more 
accurate, but physicists are accustomed to the word lattice. If our lattice does have translation 
invariance, e.g. a hypercubic lattice, we shall say so explicitly. We do, in fact, investigate translation 
invariant cases and we do discuss the cases in which the translation invariance is broken by the HF 
state. Thus, there is a special column in the tables of Sect. V for lattices that have the additional 
property of translation invariance (or some other spatial symmetry). In any case, our systems are 
always finite. 

Among the things we can prove about HF theory for a bipartite lattice is the existence of a 
phase transition from a BCS state at low temperature to a normal state at high temperature. 
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Sect. IV contains the analysis of the repulsive Hubbard model; most of the results here — but 
not all — are a transcription of the results in Sect. Ill via a particle-hole transformation. One of 
the earliest HF studies of this model was by Penn [PD]. 

One question that is peculiar to unrestricted HF theory is whether or not the orbitals (which 
are well defined even for the generalized, particle nonconserving theory) are simple products of 
spatial functions and spin functions (the latter being one of two types, either spin up | or spin 
down I). For the half- filled band and real hopping matrices we can show this to be true for both 
the attractive and repulsive models; this is one of our more complicated proofs, and it involves a 
somewhat delicate convexity argument. 

Another question concerns the uniqueness of the HF state for a finite system. Apart from pos- 
sible global gauge transformations (those describes as "broken" in Tables 1-3 in Sect. V) uniqueness 
does hold for a half-filled band as we prove in Theorems 3.12 and 4.5. 

The main thing one wishes to know about the true HF state is whether or not it is qualitatively 
correct. Prom that point of view, a main question is whether the HF state breaks the symmetries 
of the problem, and if it does so in conformity with what one believes to be the case in the 
corresponding exact quantum state. For example, the repulsive Hubbard model on a bipartite 
lattice has total spin equal to \\A\ — \B\\ (where |yl| and \B\ are the number of sites in the two 
sublattices) in a finite system ground state [LE2] (see also [LE3]) and is expected to have Neel long 
range order in three or more dimensions. This Neel state partially breaks the original translation 
invariance of the Hamiltonian (if there is any to start with) into a smaller group consisting of 
translations on each sublattice separately. This is exactly what we prove to be the case for the HF 
ground or Gibbs state. Indeed, after a suitable rotation of the spin basis, we will find the spins to 
point upward on the A-sublattice and downward on the S-sublattice This validates the predictions 
of mean field calculations for the translation invariant case in the physics literature [D JK] , [CM] , 
[FE] . As far as we know, this self-consistent antiferromagnetic A-B spin structure was assumed to 
be valid in the energy minimizing ground state, but it was never proved that this was indeed true. 
In principle, some sort of further symmetry breaking could occur. Our results show that this does 
not happen. 

Sect. V summarizes what we can prove about the breaking of symmetry in different cases. Un- 
fortunately, different combinations of conditions have to be treated separately; the basic possible 
postulates are bipartite (or non-bipartite) lattice, real (or complex) hopping matrix, repulsive or 
attractive interaction, half-filled (or not half-filled) band, translation invariance (or no invariance). 
The symmetries to be investigated are spin SU{2), pseudospin SU{2), U{1) (particle conserva- 
tion) Z2 (particle-hole symmetry) and translation invariance. For the reader's convenience, our 
conclusions are encapsulated in three tables. 

The phrase "symmetry breaking" does require a precise definition and we supply that in 
Sect. V. To us it means that the state fails to have the symmetry that the Hamiltonian has. Some 
authors, e.g., [BR] use a more restrictive definition, namely that the state also fails to belong 
to exactly one irreducible representation of the symmetry group. For example, the completely 
magnetized ground state of a Heisenberg ferromagnet breaks rotation symmetry in our sense, but 
not in the restrictive sense. As a matter of experience, however, this case is rare, i.e., usually the 
two senses agree in practice. At least this is so in the cases we can analyze completely here. In any 
event, as explained in Sect. IV, the restrictive definition requires a linear structure and, therefore, 
is not appropriate for HF theory. 

Acknowledgment: We thank J. Poelchau, H.T. Yau and A.S. Wightman for helpful discussions. 
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II. DEFINITION AND PROPERTIES OF GENERALIZED HF STATES 

Since its introduction in 1930 the notion of a Hartree-Fock ground state and positive tem- 
perature state has evolved. Our purpose in this section is to state clearly several definitions of 
these states, demonstrate their equivalence and prove some of their fundamental properties. De- 
spite years of attention to the subject it is surprising that some of the basic properties have not 
been clearly stated, much less proved. Two of these are in Theorem 2.12 (the variational principle 
[LEI]) and Theorem 2.14, which states that there are newer unfilled shells, regardless of the particle 
number. 

Il.a DEFINITIONS OF HF STATES 

The original point of view was that a HF state is a single determinant of one-particle orbitals, 
in the variables zi,Z2, ■ ■ ■ , z^, where z = {x, a) denotes a space-spin variable for one particle. The 
inner product of two functions is ((^i|i^2) = Y1I '^^2{x,cr)dx. The usual HF state is then 

a 

^UF = (iV!)-'/'Det[(^,(x,)]i<i,,<iv =: {Nl)-^/\ip, A...A^n), (2a.l) 

in which Det[(99i|(/?j)] = 1. 

This is not general enough for our purposes because we also want to allow for HF states in 
which the particle number is not conserved. After all, there are other quantities such as total spin, 
total momentum and total angular momentum that are not necessarily sharp in a HF state, and 
there is no reason why the particle number should not suffer a similar fate. In any case, when we 
go to positive temperature, the HF state should not be expected to be a pure state if it is to have 
any physical relevance. For these reasons we are going to make definitions that go beyond simple 
determinants. There is nothing new about the definitions given here but it is important for us to 
be very clear about them. 

Abstractly we begin with a one-particle Hilbcrt space Ti (finite or infinite dimensional, but 
always separable) . We define the fermionic A''-particle space to be the antisymmetric tensor product 

N times 

of H(^) := ?i a . . . a ?i for all a/" = 1, 2, 3 . . .. A simple vector in H^^^ is of the form 

/i A/2 A... A/iv := J2 (-)"-/-(i)®/-(2) ^•••®/7r(iV), (2a.2) 

permutations 

where each /j is in TL. A general vector in TL^^'^ is a (possibly infinite) linear combination of such 
simple vectors. We note that the squared norm of the above simple vector is^" 

(/i A . . . A /Ar|/i A . . . A /iv) = iV! Det[(/i|/,)]i<i,,<jv. (2a.3) 

With the identification TY^^) = H and H^^^ = C • |0) where |0) is the vacuum and (0|0) := 1 
(of course, |0) ^H),we can define the Fock Hilbert space 

jF := ^^(0) e e ^^p) ® . . . . (2o.4) 

^ Dirac notation will be used. {f\g) is the inner product of / and g, which is linear in g and conjugate linear 
in /. (/lAlgi) = {f\Ag) is the inner product of / with the vector Ag, i.e., the operator A acts to the right. It is 
important to keep this in mind when A is not self-adjoint. 
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To any vector / G H we associate a creation operator c^f) and an annihilation operator c(/), 
each acting on J^. The creation operator acts on simple vectors by 

ct(/)(/i A ... A /at) = (TV + A /i A . . . A /;v). (2a.5) 

This definition extends to by linearity, and c(/) is defined to be the adjoint operator of c^(/). 
Note that c(/)|0) = and that c(A/) = Ac(/) for any f and any A G C. (Here and elsewhere 
complex conjugation is denoted by a bar.) By this construction the creation and annihilation 
operators fulfill the canonical anticommutation relations (CAR): 

{c(/),ct(5)} := cif)c\g)+cHg)c{f) = {f\g) ■ 1, 

{c\f),c\9)} = {c{f),c{g)} = 0. (2a.6) 

Here, 1 is the identity operator on J^. We remark that the Fock space is determined by the 
vacuum |0) and a complete set of operators c^(/), c(/) that obey the CAR. Indeed, simple vectors 
can be written as 

/i A . . . A /jv = (7V!)^/2ct(/i)ct(/2) . . . ct(/^)|0). (2a.7) 

A Bogoliubov transformation of .7^ is a unitary operator W, on with the following 
special property: For each vector f E H 

d\f) := >Vct(/)>Vt = c\g) + c{h) = c\uf) + c{vf). (2a.8) 

Here g and h arc vectors in Ti and, by the linearity of W we see that there must exist a linear 
operator u on Ti. such that g = uf and there must exist an antilinear operator v on TC such that 
h = vf (antilinear means that v{Xa + b) = Xv{a) + v{b)). It is easy to check that the unitarity of W, 
which in particular implies the CAR for the operators d^f),d{f), results in certain conditions on 
u and V. Because v is antilinear these are not easy to state. One way is to choose an orthonormal 
basis /i,/2,... in Ti and define the matrix elements Vij = {fi\vfj) and Uij = {fi\ufj). In terms 
of the linear operators u and v having the same matrix elements the condition is that the linear 
operator 

acting on H (B Ti, is a unitary matrix. Here u,v denote the linear operators with the complex 
conjugate matrix elements Uij,Vij. An equivalent, basis independent statement is this: to every 
antilinear operator v we can always associate an antilinear operator such that {f\vg) = {v'^ f\g) 
for every f,ge'H. The condition becomes 

u^u + v^v = 1 = uu^ + mF 
V U + U^V = Q = UV + vu^ 

In case H is finite dimensional the converse is also true, i.e., if a pair («,v) satisfies (2a.l0) 
then there is a unique unitary operator W on satisfying (2a. 8). If Ti is infinite dimensional it 
may not be possible to find a W, even if (2a. 10) is true. This phenomenon occurs, for example, 
when -u = 0. A sufficient condition for the existence of W is given in Theorem 2.2. 

The second notion needed for generalized HF theory is a quasi-free state. In general a state, p, 
on the set of bounded operators that act on JT is a (complex- valued) linear map (i.e., for all operators 
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A and B, p{XA + B) = Xp{A) + p{B)) satisfying the conditions p(l) = 1 and piA'fA) > (which 
imphes P{A) = p{A^)). The example that will concern us most is a pure state P{A) = (^|yl-0) for 
some e Another important example is the Gibbs state p{A) = Z~^Tr[Aexp(— with 
Z = Tr[exp(— for some Hamiltonian H on T with Z < 00. 

A state p is quasi-free if all correlation functions can be computed from Wick's Theorem, 
i.e., if the operators 61,62, ■■■ , 62N are each either a c''' or a c, then p(eie2 • . • e2Ar-i) = and 

P(eie2 . . . 62n) = 2^ {-TP{6^{l)6n{2)) ■ ■ ■ P{6„{2N -l)6n{2N)) (2a.ll) 
7r 

where Yl' -r sum over permutations vr which satisfy vr(l) < vr(3) < . . . < 7r{2N — 1) and 

7r(2j — 1) < 7r(2j) for all 1 < j < A''. The right side of (2a. 11) is also known as the Pfaffian of the 
triangular array [p(ejej)]i<i<j<2Ar. In particular, we have the important formula 

^(61626364) = ^(6162)^(6364) - ^(6163)^(6264) + ^(6164)^(6263). (2o.l2) 

When this is applied later to the expectation value of the two-body potential these terms will 
correspond to the direct, the exchange and the pairing energies (see (2c. 8)). 

We remark that the quasi- free states are invariant under Bogoliubov transformations, i.e., 
if p is quasi-free and W is a Bogoliubov transformation, then the state Py\;{A) := p{yVAW^) is 
quasi-free, too. 

If p is a pure state, i.e., p{A) = {tjj\A\tp), and if the vector tp lies solely in some fixed H^^^ 
(including the possibility N = 0) then p is a quasi-free state if and only if is a normalized simple 
vector (including the possibility = |0)). This state is the usual N -particle Slater determinant 
state defined by taking expectation values with respect to the vector ipuF given in (2a. 1). 

We can define the (unbounded) particle number operator on by the formal sum 

00 

Ar= Arn(^), (2a.l3) 
Ar=o 

where II*^^^ is the projector onto the subspace H^^^ C JF. (Note that 1 = J2n^^^^-) ^ state p 
has finite particle number if 

00 

p{J\f) ■= Arp(n(^)) (2a.l4) 

JV=0 

is finite. These are the states of primary physical interest. 

A generalized Hartree-Fock state is defined to be any quasi-free state having finite particle 
number. 



Il.b ONE-PARTICLE DENSITY MATRICES 

Let p be a state and let {/i, /2, . . .} be an orthonormal basis in 7i. We define the one-particle 
density matrix (l-pdm) F to be the self-adjoint operator on'H®'H whose matrix elements are 



r 



P([ct(5i) + C(52)] [C{h,)+C\h2)\), 



(26.1) 
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where g := J^'Pkfk when g = Y1 f^kfk- Note the ordering of the operators here. The matrix 1 — F 

k k 

has a more "natural" appearance; using the CAR we find that 



r)(^^)) = P{[c{hi) + cHh,)][c\g,)+c\g,)]). (26.2) 



Note that the definition of T as an operator on 7i (B TC depends on the choice of the basis 
{/i) f2, ■ ■ ■}■ If this basis is changed then F itself changes if the antiunitary map g ^ g changes 
(which may or may not occur) . The underlying reason that F cannot be defined in a basis indepen- 
dent way is that F is not, intrinsically, a linear map on 7i©7Y; it also has an antilinear component. 
This basis dependence is the price we pay for avoiding the introduction of an abstract antiunitary 
map. The quantities we are going to compute later by means of the 1-pdm F will, however, be 
independent of the choice of the basis {/i, /2, . . .}, which we shall consider to be fixed henceforth. 

2.1 LEMMA: For any state p and any orthonormal basis {/i, /2, . . .} 

< F < 1 (26.3) 

holds as an operator onH^BH. 

Proof: Let ^ = {f, g) e H ® H he normalized. Then 1 = ||^f = ||/f + = + Hsf 
and, from (2b. 1) and the CAR, it follows that 

< {ip\T^) = p[ic\f) + c(5))(c(/) + c\g))] 

= (||/f + ||ff||>(l)-4(c(/) + ct(ff))(ct(/)+c(5))] <1. ■ (26.4) 

It is convenient to view F as a 2 x 2-matrix of operators on H. Writing 
one easily finds, using c|. ■= c^{fk), that 

{fm 

{Ua^fk) = P{cicl), (26.6) 

where the operator A is defined by 



ifmhfk) = P{cicJ, 



Note that 



iUAfk) := iUAfk). (26.7) 
= q;^ = _a (26.8) 



(where := at). 

In view of (2a. 11) and the density of polynomials in the algebra of observables, the 1-pdm F 
of a quasi-free state p uniquely determines p. More importantly, any F is the 1-pdm of a quasi-free 
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state /3, as we shall show in Theorem 2.3 below. We shall, however, restrict attention to finite 
particle states. It is easily seen that N = Ylk ^k^k and that, for any state, 

k 

First, we give the relationship between the 1-pdm F of a state p and the 1-pdm Fyy of the 
transformed state Pw, Pw{A) := p(WAW^), assuming W is a Bogoliubov transformation. Using 
(2a.9) we can write the Bogoliubov transformation (2a.8) as 



dk\ sr^ I ^ki ^fci \ ^ Ci\ x-^,„.+ . /c. 




Note, as indicated above, that each {W^)ki is a 2 x 2 matrix. We then find that 

{i-r^)ki=('^'^f^ '^"^f^] =5:(Tyt),,(i-r),,(w^),,. (26.10) 



,P{dldi) Pididi), 



Thus 

Tw = W^TW. (26.11) 

We now give a sufficient condition for the operator W in (2a.9) to represent a Bogoliubov 

transformation in JF. 

2.2 THEOREM: Any unitary operator W of the form (2a. 9) that satisfies the condition 
Tv[vv'^] < 00 always corresponds to a Bogoliubov transformation W on with 

S{f) := Wct(/)Wt = c^uf) + c{vj). (26.12) 

Moreover, W|0) = |^) is the state 

1^) = n{(l - CO'/' + Cy"c\h,)c\k,)}f[cH9im. (26.13) 

i i=l 

Here, (1,(^2, denote the eigenvalues of vv^ in the interval (0,1) counted with half their mul- 
tiplicity (these eigenvalues are all evenly degenerate), and r is the multiplicity of 1 as an eigen- 
value of vv'^ . The vectors 91,92, ■■■ ,9r, hi, ki,h2,k2, ■■ ■ form an orthonormal family of eigenvec- 
tors of vv^ with gi, g2, ■ ■ . , gr being the eigenvectors of eigenvalue 1. The pair of vectors hi, and 
ki := (Ci ~ C,f )~'^^'^'uv'^ hi are eigenvectors of vv^ of eigenvalue Q. 

Remark: We refer the reader to [AH] where the first statement of the lemma, together with 
its converse (which is not needed here), is proved ([AH], Theorem 7). However, we are not aware 
that the explicit formula (2.22), which we do require, is readily accessible. 

Proof: The unitarity of the operator W implies the following conditions similar to (2a. 10) for 
the linear operators u and v: 

u^u + v^v = 1 = uu"^ + vv^ 

+ _+ _+ (2&-14) 

v^u + wv = = uv' + vw . 
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Thus, < vv^ < 1 as an operator on Ti with purely discrete spectrum, thanks to Tr[i;i7^] < 00. 
Furthermore, if /i is a normahzed eigenvector of vv^ with eigenvalue < ^ < 1 we find, using 
(2b. 14), that A; = (C - O'^^^uv^h satisfies 

vv^k = (c - c^y'^^^vv^uv^h = -(c - c'^y^^^vu^vv^h 

= -c(c - c^um = ck (26.15) 

and thus k is also an eigenvector of vv^ . Moreover, 

{h\k) = ic- er'^Hh I uvm) = -(c - e){h I vu^h) 

= -(C - C^){uv''h \h) = -{h\k) (26.16) 

and hence {h \ k) = 0. Likewise, wc sec that {k \ k) = 1. Iterating the map from h to k will not 
produce more eigenvectors since {(^ — C^'^)~^/^uv^ k = {( — (^)~^uv^uv'^ h = {( — (^)~^uu^vv^h = —h. 

We can thus find an orthonormal basis for H of the form pi, . . . , ^ij ^1, ^2, ^2, • • • , 
/i , /2 ! • • • ) where gi,g2, ■ ■ ■ and li,l2, ■ ■ ■ are the eigenvectors of vv^ of eigenvalue 1 and 0, respectively. 
Another orthonormal basis is given by v^^i, ■ ■ ■ jV^dn Ci ^^'^v^hi, ^^'^v^ki, . . . , 
u^li,u^J2, ■ ■ ■■ To prove this we first note, as is easily seen, that this is an orthonormal fam- 
ily. We then note that if / G H is orthogonal to all members of this family then v^vf = 
(because vf is orthogonal to Qi, . . . ,gr,hi,ki, . . . but not necessarily to Zi, . . . ; however anni- 
hilates 7i, . . .). Hence v'^uf = —v)vf = 0. Thus uf is an eigenvector for vv"^ of eigenvalue 0, i.e., 
uf G span{/i, Z2, • • •}, but then / = (1 — v^v)f = u^uf G span{uUi,uH2, ■ ■ .}, and thus / = 0. 

Using this latter basis we define, in agreement with (2b. 12), 

d{v^gi) := c{uv%) + c\vv^gi) = c\gi) i = 1, . . . ,r, (26.17) 

where we have used uv^'g^ = —vu^g^ = 0, which follows from uu^gi = (1 — vv^)gi = 0. We also 
make the definitions, for i = 1, 2, . . ., 

d{C-'/h%) := C-'^^c{uv%) + Q^'^^vv'^h,) = (1 - Cif'^c{h) + Cy'ct(/i,), (26.18) 
^^(Cr'^Vfe,) := C-^'\{uv%) + C-''^c\vv^h) = -(1 - Cif'\{hi) + Ci'^c\h), (26.19) 
d{v)li) := c{uuHi) + c\vu%) = c(li). (26.20) 

We shall now show that (2b. 13) defines a normalized vector l'^) in .7^ annihilated by all the 
operators in (2b. 18-20). To show that the somewhat formal expression on the right side in (2b. 13) 

defines a vector in we expand the (possibly infinite) product, thereby arriving at a (possibly 
infinite) sum of (possibly infinite) products. Each term in this sum that contains a product of 
infinitely many c^'s is zero since it will also contain a product of infinitely many Cl^^ ^ ^-iid Q ^ 
as i 00. Hence, the sum is at most a countable sum of orthogonal simple vectors in J^. An 
appropriate truncation of this sum gives 

iV'iv) := U(l-Ci)'^'Ul(^-Ci)'^' + Cl^'cHhi)cHk,)]flcH9m- (2b.22) 

i>N i<N i=l 
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Wc note that {iIj^ \ iPn) = YliyN^^ ^ Ci) ^ 1 ~ X]i>Af d- which is non-zero and converges to 1 as 
N —>■ 00 since Yli Ci < 00. li M > N then \tpN) is orthogonal to \tpM — i^n) and 

(V'M - V'iv I V'M - V'iv) = (V'M I V'm) - (V'iv I V'iv) = Yli^- Ci) ~Yl{^-Ci)^0 (26.23) 

i>M i>N 

as M ^ 00. Hence there is a normalized G ^ such that iV'iv) — * \ip) as ^ 00. It is easy to 
check that any operator in the list (2b. 18-20) will annihilate IV'at) if N is large enough. Thus IV') is 
annihilated by all the operators in (2b. 18-20), and hence by any operator d{f) := c^{uf) + c{vf). 
Denoting d^. := c\ufk) + c{vfk) (recall fj^ = fk) we define W : T ^ T hy 

Wc\...ci\0)=4...di\i;) (26.24) 

for all n > 0. The unitarity of W impHes that satisfy the CAR and hence that W is an isometry. 
It is also clear that VV satisfies (2b. 12). 

To show that W is unitary we first observe that |0) can be written as a linear combination of 
vectors of the form (2b. 24). This is easily seen by interchanging the roles of |0) and and of W 
and W"^ in the argument which led to the construction, (2b. 13), of Prom this it follows that 
all simple vectors c\ . . . cJ^jO) can be written as combinations of the vectors of the form (2b. 24). 
Indeed, from (2b.l4), 4 = ct((nut + vv'i)fk) + c{{uv'' + TO^)7fe) = (i^ut/fe) + div^fj^). Thus W is 
invertible and hence unitary. ■ 

2.3 THEOREM: Let < F < 1 be an operator on H ® H of the form (2b. 5), subject to 
(2b. 8), and assume furthermore that Tic[y] < 00. Then there exists a unique quasi-free state p with 
finite particle number such that T is the 1-pdm of p. 

Remark: Suppose we are given a 7 satisfying < 7 < 1. We can setr=^Q 1*^7^' 

set a = and, according to Theorem 2.3, we can extend F to a quasi-free state p. In other words, 
given 7 we can find a particle conserving quasi-free state having this 7 as its 1-pdm (whether or 
not Tr7 is an integer). By a particle conserving state we mean a state (such as a Gibbs state for 
a particle conserving Hamiltonian) that is a convex combination of states, each having a definite 
particle number, i.e. there are no cc or c^c^ matrix elements in this state, which is the same thing 
as saying a = 0. However, unless 7 is a projection (i.e., a usual HF state) F cannot have a definite 
particle number, i.e. p(A^^) > p{MY . In fact, defining p{J\f'^) = Y,n -^^P(7r^^^) = Y^k i Pi^^k^kcjci), 
we can use (2a. 12) to compute 

PiAf) = PiAff + J2 {Vki + ^Ik^ik) +J2^kk = Pi^f + ^[7 - 7'] + ^[«^«], (26.25) 

from which we see that p{J\f^) = p{MY requires both a = and 7^ = 7. 

Proof: Since Tr[7] is finite it is clear from the form in (2b. 5) that although Tr[F] may be infinite, 
which is the case when dimH = 00, we have Tr[F(l — F)] < 00. Thus, there is an orthonormal 
basis of eigenvectors for F(l — F). If ip is an eigenvector for F(l — F) of eigenvalue p then so is Tip 
and, since T'^(p = T(p — ptp, it follows that F leaves invariant the subspace space {ip,rip}, which is 
at most 2-dimensional. We conclude that there is an orthonormal basis of eigenvectors for F. 

If (y? = / © (7 is a normalized eigenvector for F of eigenvalue A then, by (2b. 5) and (2b. 8), we 
find that ip = g ® f is a normalized eigenvector for F of eigenvalue (1 — A). Thus, we can find a 
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unitary W onH(B'H oi the form (2a. 9) sTich that (using the basis /i, /2, • • • for both copies of H) 
the four blocks of the transformed F have the form 



/Ar 



w^rw 



\ 



l-Ai 



1- Ao 



V 



(26.26) 



•••/ 



with < Ai < 1/2 for i = 1, 2, . . .. Using this and Tr[r(l - T)] < 00, we obtain J^i Ai < 00. 

We shall now prove that W satisfies the condition of Theorem 2.2, i.e., that Tr[t;u^] < 00 and 
hence that W corresponds to a Bogoliubov transformation W on T. Indeed, from (2b. 26) we know 
that the upper left block of the matrix VF^TW has finite trace, i.e.. 



Tr[tt^7ti + v^a^u + u^av + v^l — 'y)v] < 00. 



(26.27) 



Since < F < 1 wc have that F^ < F. The upper left block of this inequality reads 7^ + aa^ < 7. 
Thus Tr[aa'''] < Tr[7] < 00. By Cauchy-Schwarz we estimate 

Tr[ti'''7'U + v^a^u+u^av + v^(l — j)v] — Tr\y^v] 

= Tr['u^7M + v^a^u + u^av — v^jv] 

> -Tr[u'U^7] - 2Tr[v^v]^/'^Tr[uu^aa^]^/'^ - Tr[vv^j] 

> -2Tr[7] - 2Tr[w%]^/^Tr[7]^/^ (26.28) 

We conclude from (2.31) that Ti[v''v]-2Tr[v'<v]^/'^Tr[j]^/^ < 00 and hence that Ti[vv^ = Ti[v^v] < 
00. 

If we prove that the diagonal matrix W'^'FW is the 1-pdm of a quasi-free state p then we know 
that F is the 1-pdm of the state p with p{A) = p{W^AW). We may therefore assume that F is 
itself diagonal of the form (2.26). 

Let Ho be the projection onto the subspace of JT on which Yli-\i=o ~ ^■ 

For each i such that Aj > choose Cj to satisfy 



(1 + exp(ej)) ^ = Aj, 



(26.29) 



(note that < ej < 00, since < Aj < 1/2) and consider the following (possibly unbounded) 
operator H on JT. 

4ci. (26.30) 



i:Ai#0 



We shall now prove that the operator 

G := Uo exp(-jy) 

has finite trace on and that the state 

P(A) = Tr[G]-^Ti[AG] 
13 
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is quasi-free and has F as its l-pdm. It is easy to see that the trace of G is 



Tr[G] = n + exp(-ei)] = - A^)"^ < 00, (26.33) 



since ^ Aj < 00. 



The operator G looks peculiar but, by introducing the operator H' = Yl A^i which com- 

i:Ai=0 

mutes with iJ, we can write the state p as a limit of Gibbs states: 

P{A)= lim Zr^Tr[Aexp(-/?if'-if)]. (26.34) 

It is well-known (see [GM] for a simple proof) that the Gibbs state for an operator of the form 
X] Cjcjcj is quasi-free. Hence we conclude (by taking the limit /3 00) that p is quasi-free. 

i 

The fact that V is the l-pdm of p follows from the computation 

^((4 + c,)(c„ + 4)) = Tr[G]-iTr[(4 + c,){c^ + 4)G] 

= 4m1V[G]-^Tr[4c;,G] +5;„TV[G]-^Tr[cjc|G] 

= 11(1 - Ai) I 5km\k{l - \k)-^ n (1 - ^^)"' + n (1 - ^i)"' 

= Skm>^k + Sinil-Xi), (26.35) 

where we have used the fact that exp[— Cfc] = Afe(l — A^)"^ if A^ ^ 0. Finally p has finite particle 
number since P{J\f) = Tr[7]. 

The uniqueness of p follows as in the remark after (2b. 8): The l-pdm of a quasi-free state 
determines the state. ■ 

We call an operator F admissible if it satisfies the properties in Theorem 2.3, i.e., is of the 
form (2b. 5) subject to (2b. 8) with Tr[7] < 00 and < F < 1. We then have that: 

F is admissible if and only if it is the l-pdm of a generalized Hartree-Fock state. 

In the above proof we not only proved the existence of a quasi-free state having F as its l-pdm, 
we also gave the explicit form of p. To make this more explicit we introduce the following notion: 

Quadratic Hamiltonian: A self-adjoint operator H (hounded or unbounded) on T is said 
to be a quadratic Hamiltonian if the unitary operators W(t) := exp{iHt) are Bogoliubov transfor- 
mations for all t. 

If is a quadratic Hamiltonian there correspond operators W{t) on H©7Y of the form (2a.9) 
corresponding to the Bogoliubov transformations W(t). Since the anticommutator satisfies 

{c\h2) + c{h,), eMiHt) (ct(5i) + cig,)) exp(-iFi)} = W{t) (^^'^ ^ 1, (26.36) 

it follows that W{t) is a strongly continuous one-parameter group of unitaries on TC (B 7i. Hence 
there is a self-adjoint operator A (bounded or unbounded) such that W{t) = exp(zAi). 
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The operator A has the block structure 



, a 6 , 

A=[,, _ , (26.37) 



6t 



-a 



where a and b arc operators on Tl with = a and = —b. 

We call^ the operator A the first quantization of H and we cah H a second quantization 
of A. Notice that from (2b. 36) A is determined uniquely by H. If i? is a bounded operator on T 
then A is bounded on'H®'H and by differentiating (2b. 36) we can, in this case, write 



{c\h2) + c{hr),[H,c\gr) + c{g^)]} = 



A[^^])l. (26.38) 



Here wc have introduced the commutator [Ki,K2\ '■= K1K2 — K2K1. 

The operator H , however, is only determined by A up to addition of a multiple of the identity. 
If A is bounded we may write the unique second quantization, H, of A satisfying (0|ff|0) = in 
terms of the matrix elements of a and 6, as 

If is a quadratic Hamiltonian and if W is a Bogoliubov transformation then WHW"^ is also 
a quadratic Hamiltonian. If is a second quantization of A then WHW^ is a second quantization 
of WAW^ , where W is the unitary given in (2a.9). 

The proof of Theorem 2.3 implies the following result about the structure of quasi-free states. 

2.4 LEMMA: Let p be a quasi-free state with finite particle number, i.e., p{Af) < og (in terms 
of its 1-pdm this means Ti['y] < 00). Then there exist two commuting quadratic Hamiltonians H 
and H' (possibly H = or H' = 0, but not both) such that 

p{B) = lim Tr[exp(-/?JT' - jy)]-^Tr[S exp(-/?jy' - H)]. (26.40) 



This means that the state p is a product of the ground state (zero temperature state) for H' 
and the Gibbs state for H. Thus if iT' = 0, p is a Gibbs state and if = 0, p is a pure state. In 
the next two lemmas we discuss quasi-free Gibbs states and quasi-free pure states in more detail. 

2.5 LEMMA: If A is an operator onU^U of the form (2b. 37) then 

r := (1 + exp(^))-i (26.41) 

is of the form (2b. 5) subject to (2b. 8). Furthermore, if this T satisfies Tr[7] < 00, so that it defines 
a quasi- free state p by Theorem 2.3, then this p is given by 

P{B) = Tr[exp(-J?)]-^Tr[Bexp(-iy)], (26.42) 

where H is any second quantization of A. 

^ Our terminology is far from being conventional, but it is descriptive. 
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Conversely, if p in (2b. 42) is a quasi-free state of finite particle number with 1-pdm T, then 
H is quadratic with first quantization satisfying (2b. 41)- 
Proof: Let U be the unitary 

U=[l J) (26.43) 

on neH. Then UAU^ = -A, since A has the form (2b.37). Hence, UFW = 1-T which proves 
that r is of the form (2b.5) subject to (2b.8). 

If r satisfies Tr[7] < 00 it follows from the proof of Theorem 2.3 that there is a Bogohubov 
transformation W with corresponding W such that the VF^Fl^^, which is the 1-pdm of the trans- 
formed state Pw, is diagonal in our chosen basis. We denote by Ai, A2, . . ., those eigenvalues of 
W^TW in the interval [0, 1/2) together with half the eigenvalues equal to 1/2 (if any). The other 
half of the eigenvalues of VF^rVF arc then given by 1 - Ai, 1 - A2, . . .. The operator W^AW 
is also diagonal with the first half of the eigenvalues given by 61,62, .. . and the second half by 
—61, —62, . . ., according to the definition (2b. 41). 

Since 

exp^i^6fec|.Cfc^c|'exp^-z^6fcc|.Cfe^ = exp{iei)cj (26.44) 
k k 

we see that the second quantizations of W^^W are of the form Ht = X^^fcC^Cfe + rl, where r is 

k _ 

any real number. Since the Gibbs states are independent of r we see that all the operators H^- 
define the same Gibbs state as the operator in (2.33), i.e., the state Pyy. (As in (2b.42) we are 

referring to Gibbs states with the inverse temperature (3=1.) 
Since Pw is the Gibbs state for Ht, the state p is given by 

p(B) = PyviW^BW) = Tr[exp(-#^)]-iTr[W^SWexp(-#^)], (26.45) 

which agrees with (2b.42) if if = WHt-W^. Such an H is however, a second quantization of 

WW'fAWW'f = A. 

The converse statement is also a simple consequence of the proof of Theorem 2.3. One only 
has to realize that the state p in (2b. 42) uniquely determines the operator H. ■ 

2.6 THEOREM: A quasi free state p with finite particle number is a pure state P{B) = 
{ip\B\ip) if and only if the 1-pdm T with Tr^y] < 00 is a projection onTi^Ti, i.e., = F. 

In terms of F the vector is of the form (2b. 13) but this time with gi, . . . ,gr,hi,ki, . . . 
being orthonormal eigenvectors ofj; the vectors gi,...gr with eigenvalue 1 and the pair hi, ki := 
— {Ci — Ci)~^^'^CKhi with eigenvalue d, where < ("i < 1- 

Proof: Let p be a pure state. Since p is uniquely determined by its 1-pdm F, we may assume 
it to be of the form (2b. 32). The purity of p, then, is equivalent to G being rank one and hence 
Aj = for all i, which is equivalent to F to be a projection. 

Again from the proof of Theorem 2.3 the diagonal form of F is 

W^FW=(^^ J^. (26.46) 



Hence 



?Wt=fr! (26.47) 
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and we have 7 = vv^ and 



-uv 



T7t 



. Since W^TW is the 1-pdm of the pure state 



corresponding to the vacuum, T is the 1-pdm of the pure state corresponding to \iJj) = W\0). Here 
W is the Bogohubov transform defined by W. The last statement of the lemma now follows by 
comparison with Theorem 2.2. ■ 

Using this lemma we can find a basis for 'H®'H where the blocks 7 and a of F take a particularly 
simple form when F is a projection. In fact, if we choose the basis consisting of Qi 0, Zj 0, hi ki 
and 7j, fcj hi for aU i = 1, 2, . . ., we find (with = (Ci - Q)'^^'^) 



I 



\ 



\ 



■J 



and a 



I 



\ 



CKj 
-OLi 



(26.48) 



2.7 LEMMA: Let p he a pure and quasi-free state of finite particle number with 
1-pdm r. Then p{AP) < 00 and it is given by (2b. 49) with a given in (2b. 5). 

P{^f^)- p{Mf = 2TT[a^a]. 



(26.49) 



This equation shows that p is not necessarily a fixed particle number state. 

Proof: Since p is pure, F is a projection and hence 7 = 7^ + ota^ which together with (2b. 25) 
implies (2b.49).B 

From Lemma 2.7 we see that a generalized HF state with 1-pdm F has conserved particle 
number if and only if the component a of F vanishes. We call a generalized HF state for which 
this holds a normal state. 



II.c THE GENERALIZED HF FUNCTIONAL 

In this section we shall introduce the generalized Hartree-Fock approximation for a self-adjoint 
operator H on Fock space of the form 

H = h-i-V, (2c.l) 



where /i is a quadratic (particle number preserving) operator 

h = ^hijclcj, 



(2c.2) 



and y is a quartic operator (again particle number preserving) 

^ — \ 5^ ^kl;mnc\.c\cnC^. 
k,l,m,n 



(2c.3) 



By (2c. 2) the operator h is defined in terms of matrix elements hij of a self-adjoint operator 
h on H. The operator h is the restriction of h to the one-body space H^^^ = H, i.e., 



h^j = {fi\h\fj) = {0\cM\0). 



(2c.4) 



By (2c.3) the operator V is defined in terms of matrix elements Vki-mn of a self-adjoint operator 
V on H'SiH. Note that we are not restricting V to be an operator on the antisymmetric two-body 
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space H*^^) =HA'H. The restriction of V to the antisymmetric subspace is equal to the restriction 

{0\ciCkVcl^cl\0) = \{Vkl-mn + Vlk-nm " Vkl-nm " Vifc;mn) 

= \{fkAfl\V\fmAfn). (2c.5) 

The operators V and h defining V and h may be bounded or unbounded. We shall, however, 
assume that the operator H is bounded below This is the case if, for example, h and V are bounded 
below. One way to ensure this is to assume that h and V are bounded below by (negative) operators 
of finite trace (it is not enough to assume that h and V are bounded below in order to have h and 
V bounded below.). 

The expressions (2c. 2-3) are somewhat formal. A more precise definition can be given as 
follows. On each T^^-'^^ we can, in the obvious way, define the sum h^^^ = X^j-^ hi of A'' commuting 
copies of h and the corresponding sum V'^^^ = Ylii<i<j<N ^ commuting copies of V. Then 

In discussing Hamiltonians of the form (2c. 1) we have two particular examples in mind. The 

first is the Hubbard Hamiltonian defined in Sect. 111. The second is the atomic Hamiltonian with 
TC being the square integrable (spinor valued) functions on R"^ and h = —(fi^/2m)A — Z/\x\ — /x, 
where /x < is a chemical potential and V = e^\x — y\~^ . Both h and V are independent of 
(diagonal in) spin. The Hubbard Hamiltonian is a bounded operator (in fact a finite dimensional 
matrix) while the atomic Hamiltonian is unbounded but bounded below. [If < then h and V 
are bounded below by operators of finite trace, but if /i = 0, h will not be bounded below by such 
an operator (e.g., the negative eigenvalues of hydrogen are not summable). This, however, is not 
a real problem because the operator H is still bounded below.] 

For both the Hubbard Hamiltonian and the atomic Hamiltonian we shall be interested in the 
ground state and its energy. The ground state is simply the (maybe not unique) state Po (with 
finite particle number) for which P(i^) takes on the smallest possible value - the ground state 
energy - provided this smallest value is attained for some state at all. Otherwise the ground state 
does not exist. 

If H is unbounded the expectation P(-ff) is not necessarily well-defined. If, however, H is 
bounded below we can define P{H). This is easy to see for states that can be written as P{B) = 
Ti[GB] for some positive operator G of finite trace on Fock space (this is not true for all states, but 
we are only interested in states for which it holds). If H is bounded from below, we can without 
loss of generality assume that H is positive. Then Tr[GH] is, when expanded in the eigenvector 
basis for G, an infinite sum of positive terms. This sum then defines p{H) (possibly to be +00). 

The expected number of particles in the ground state is p{M). Since both Hamiltonians 
are particle number preserving there is a ground state p with a fixed number of particles, i.e., 
p{M'^) = p{M)'^. The number of particles TV := p{M) is thus an integer and p is, in fact, a state 
on H^^) (i.e., p(HW) = 1). 

Alternatively to specifying a chemical potential p, we could also have specified the number N 
and then considered the problem on The equivalence of the two descriptions by Legendre 

transform (equivalence of the grand canonical and canonical ensembles) requires that the ground 
state energy is a convex function of the particle number. While this is believed to be the case there 
is, to the best of our knowledge, no rigorous proof of this fact in the two models discussed. In the 
grand canonical picture the ground state energy is a concave function of the chemical potential, 
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but as long as we do not know the convexity of the canonical energy as a function of wc cannot 
assert that the two energy functions are Legendre transforms of each other. Here we shall mostly 
work in the grand canonical framework, i.e., specify a chemical potential, except at the end of the 
section where we shall discuss the canonical picture when V is assumed to be positive. 

In addition to the ground state energy we shall also be interested in the grand canonical Gibbs 
states p{B) = Z~^Tr[B exp(—l3H)]. The Gibbs state, however, is not well defined for the atomic 
Hamiltonian since the operator exp(— /3i7) will not have finite trace in this case. 

The object of study in this section is not the real ground states and Gibbs states but rather 
their (generalized) Hartree-Fock approximations which we shall now define. 

The Hartree-Fock approximation to the ground state is simply the generalized Hartree-Fock 
state with least possible energy. By Theorem 2.3 there is a one-to-one correspondence between a 
generalized HF state p and its 1-pdm T. We may therefore define the generalized Hartree-Fock 
energy functional, 

e{T)=p{H), (2c.6) 

on the set of all admissible density matrices. 

The energy of a generalized Hartree-Fock state can be computed in terms of the 
1-pdm r as follows. The expectation of the quadratic part is p{h) = Tr[h'y]. In computing the 
expectation of the quartic part we apply (2.12) and obtain 

Pi^) = ^ Yl ^kl;mn{jmk^nl-^mnnk+<^lk<^mn) ■ (2c.7) 
k,l,m,n 

The operator G^^) on H H with matrix elements gJ^^^ = 7,„fc7„/ - ImHnk + «L«mn ^as finite 
trace. In fact, if we choose V = I in (2c.7) we obtain V = J\f{M-l) and Tr[G(2)] = p{J\f^) - p{J\f) 
which is finite by Lemma 2.7. 

Equation (2c.7) states that p(F) = ^Ti:[VG^% Both Tr[VG(2)] g^^^j TT[h^] are well-defined 
since V and h are bounded from below and G^^^ and 7 are positive operators of finite trace. 

As mentioned after (2a. 12) the three terms in (2b. 9) are called respectivley the direct energy, 
the exchange energy and the pairing energy. 

We can thus write the HF energy functional as 

f (F) = Tr[/l7] + 5 ^ Vkl-mn (imklnl - Imnnk + "L'^mn) • (2c.8) 

k,l,m,n 

The Hartree-Fock energy is given by 

E^^ := inf {^^(r) | F is an admissible density matrix} . (2c.9) 

As we have proved the one-to-one correspondence between quasi-free states and admissible 
density matrices, (2c. 9) is evidently equivalent to 

E^^ :=mi{p{H) \ p is a quasi-free state} . (2c. 10) 

We shall not discuss, in general, whether the infimum in (2c. 9) is attained. For the Hubbard 
Hamiltonian, however, it is clearly the case that the infimum in (2c. 9) is attained since the set 
of admissible denisity matrices is a compact subset of a finite dimensional space. In case of the 
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atomic Hamiltonian it is also true that the infimum is attained. This result was proved in [LS], 
where it was assumed that a = 0, but this follows from Theorem 2.11 below). A 1-pdm for which 
the infimum (2c.9) is attained defines a HF ground state. 

To define the finite temperature HF Gibbs state we must introduce the entropy of a 
quasi-free state: 

S{T) := -iTr[rinr] - iTr[(l - T) ln(l - F)] = -Tr[rinr]. (2c.l0) 

The last equality in (2c. 10) holds because F and 1 — F are unitarily equivalent (cf. proof of 
Lemma 2.5) and F has real eigenvalues. Notice that by Theorem 2.6, -S'(F) = if and only if F is 
the 1-pdm of a pure state. 

We define the Hartree-Fock pressure functional Vjs at inverse temperature /3 as 

-rfs{r) = e{T)-r's{T). (2c.ii) 

The Hartree-Fock pressure is defined by 

V^^iP) = sup{7-*/3(F) I F is an admissible density matrix} . (2c.l2) 

As explained above we only consider positive temperature in case of the Hubbard Hamiltonian 
and as for the energy it is then clear that the supremum is attained. A HF Gibbs state is defined 
by a 1-pdm. maximizing (2c. 11). 

In the next lemma we show that the HF energy gives an upper bound to the true energy and 
the HF pressure gives a lower bound to the true pressure. 

2.8 THEOREM: We have the inequality 

^HF > ._ infp^H), (2c.l3) 

where the infimum is over all states p (not just HF states). If h and V are bounded 

V^^iP) < := 0-'^ lnTr[exp(-^JT)]. (2c.l4) 



Proof: The inequality (2c. 13) is obvious since E^^ is defined according to (2c. 10) as an infimum 
over a restricted class of states, namely the generalized HF states. 

Inequality (2c. 14) is more complicated. The aim is to show that for any generalized HF state 
p with 1-pdm F we have 

exp{-f3p{H) + S{T)) < Tr[exp(-/3iJ)]. (2c.l5) 

According to Lemma 2.4 any generalized HF state is a limit (in the sense of convergence of 
expectation values of bounded operators) of quasi-free Gibbs states. Moreover, it is clear that the 
entropies of the approximating Gibbs states are greater than the entropy of the limiting state. We 
can therefore assume that the quasi-free state p in (2c. 15) is a Gibbs state. 

Since p is a Gibbs state we can define an operator ^ on ?i © 7^ as in (2b.41). This operator 
is then the first quantization of a quadratic operator Ha and 

P{B) = Tr[exp(-/3/tA)]"^Tr[Sexp(-/3/tA)]- (2c.l6) 
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To specify Ha uniquely we assume that (0|/?.a|0) = 0. 
By the Peierls-BogoUubov inequahty [TW] we infer 

TV[exp(-/?if)] =Tr[exp(-/3(/iA + H- Ha)] 

>Tr[exp(-/3/iA)] exp(-/3p(/f - /i^)). (2c.l7) 

The inequahty (2c. 14) fohows if we can prove 

5(r) = (3p{hA) + In (Tr[exp(-/?/i^)]) . (2c.l8) 

By diagonahzing the operator A as in the proof of Lemma 2.5 we find 

Tv[eM-PhA)] = + exp(-/?efc)), (2c.l9) 
k 

where 61,62, ■■ ■ again denote the positive eigenvalues of A and half the zero eigenvalues (if any). 
We also find 

p{hA) = Tr[exp{-phA)]-^Tr[hAeM-PhA)] = e^fl + exp(/3efc)]-^ (2c.20) 

k 

Using (2b.41), (2c.l9) and (2c.20) we finally obtain that the right side of (2c.l8) is 

5^ In (1 + exp(/3efe)) + ^ /3efe [ (l + exp(/3efe)) ~' -l]=^Ti[- ln(r) - (1 - T) ln{r-' - 1)] , 
k k ^ ' 

(2c.21) 

which is exactly ^(r).^ 

It is always possible to choose the quantum ground state to be a pure state. The same is true 
in generalized HF theory. This follows from the next lemma when compared with Theorem 2.6. 

2.9 THEOREM: The infimum of £ over all admissible density matrices agrees with the 
infimum over all admissible projections, i.e., 

E^^ = inf{£"(r) I r is admissible and = T}. (2c.22) 

Proof: We shall show that for any admissible F there is a projection Tq such that ^^(ro) < ^(r). 
For this purpose it suffices to approximate h and V by bounded operators — for which there is 
obviously no difficulty with convergence of the following sums. 

For any admissible T, as in the proof of Theorem 2.3, we can find a Bogoliubov transformation 
W with corresponding W such that VF^FVF is diagonal. If T is the l-pdm of the state p then ly^riy 
is the l-pdm of the transformed state Pw. We have the relation ^^(r) = P{H) = Py\;{W^ HW). 
The transformed Hamiltonian W^HW is also a sum of a quadratic operator W^hW and a quartic 
operator W^VW, but these are not necessarily number preserving. If we (anti)commute all to 
the left (normal ordering) we obtain 

W^HW = Y^hA^J + I E ^kl;n.n44^nCm + ^1 + R, (2c.23) 
i,j kl;mn 
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where h-j and VJ.^.^^ are new matrix elements and k is some constant. The operator R contains 
all particle non-conserving terms of the form 

cc, c^c^, cccc, c^ccc, c^c^c^c, c^c^c^c^. (2c.24) 

Since Pw{R) = 0, we obtain 

£{r) = Pw(>V^i?>V) = K + Y^hiiXi + ^ - Vki.ik)Xkk, (2C.25) 

i kl 

where Ai, A2, • • • , are as usual the eigenvalues of F smaller than 1/2. 

The important fact to observe about the expression in (2c. 25) is that, although it is a quadratic 
form in Ai, A2, . . ., it is linear in each variable. Hence d£{T)/dXj is independent of Xj. Therefore, 
we do not increase the energy expectation by replacing Xj\ipj){ipj\ + (1 — Xj)\U'ipj){U'(pj\ in T by 
\<fj){^j\ in case d£{T)/dXj < and by \UTpj){UTpj\ otherwise. Proceeding this way, we arrive at a 
1-pdm To with energy no greater than before and with all eigenvalues either or 1. This means 
that To is a projection. ■ 

If P is a generalized HF ground state for H we define a corresponding HF mean field 
Hamiltonian. It is the following quadratic Hamiltonian written in terms of the blocks 7 and a 
of the l-pdm F of p. 



i,j kl;mn 



lkmA^n~^ lnl^\^m lml^\^n Ink'^l^m 



+ ^Ik^n^m + (^mn'^k^l 



(2C.26) 



A HF ground state of H is self-consistent in the sense given in the next lemma. 

2.10 LEMMA: If p is a HF ground state for the Hamiltonian H, i.e., a HF minimizer, then 
p is a true (not just HF) ground state for the Hamiltonian Hp. 

Proof: We must show that p{Hp) < p'{Hp) for any state p' with finite particle number. 

From the l-pdm V of p' and the 1-pdm F of p we can for < t < 1 construct a new 1-pdm 
Ft = (1 — t)T + tV . Then F^ is admissible and since F is a minimizer for £ we have 

< = p'{Hp) - piHp), (2C.27) 

dt t=o 

which proves the claim. It is important here that since Hp is quadratic, P'{Hp) depends only on 
the 1-pdm of p'. ■ 

Although, we can always find a state with fixed particle number among the quantum ground 
states for the Hamiltonian H, this may not be the case for the generalized HF ground states, as 
discussed in the introduction (and proved for the attractive Hubbard Hamiltonian in Sect. III). 
Here there need not be a normal ground state. 

It is often stated in the literature that if the HF ground state is not a normal state, (i.e., it is 
a BCS state) the potential V must have a negative component. It is now very easy to state this 
precisely. 
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2.11 THEOREM: // the operator V is positive (semi) definite onH^H 

E^^ = mi{£{y) I r is admissible and normal, i.e., a = 0}. (2c.28) 

Moreover, ifV is strictly positive (i.e., positive definite) then any ground state (if it exists) must 
be a normal state. Likewise, we have for the pressure 

V^^iP) = sup{'P^(r) I r is admissible and normal }. (2c.29) 

Proof: That V is positive means that for all g & H 'Si H we have 

X] Vkl;mngkl9mn > 0. (2c.30) 
feJ;mn 

Strict positivity means that if 5 7^ then (2.67) is a strict inequality. The pairing energy is exactly 
of the form (2.67) hence 

Yl ^kl;mn(^lk(^mn = Yl ^kl;mn^kiamn > 0. (2c.31) 

kl;mn kl;mn 

Assume F is admissible with non-zero a and form a new operator F by replacing a by zero. 
It is clear that F is still admissible and by (2c. 31) that £{T) is (strictly) smaller than £{T) if V is 
(strictly) positive. To prove the result on the pressure we notice that the 1-pdm T' obtained by 
changing a to —a is unitarily equivalent to F. Hence, since the entropy is a concave function we 
find 

SiV) = S(|F + IF') > i5(F) + i5(F') = S{r). (2C.32) 

■ 

In the remainder of this section wc shall specialize to the case of positive V and discuss the 
canonical picture. Thus, instead of a chemical potential we now fix the value Tr[7] = N. Since we 
have a = it is enough to consider the component 7 of F. 

The discussion is of special interest in the atomic case where the potential V is strictly positive. 

There is a version of Theorem 2.9 for the canonical case of fixed particle number A'^. We state 
it below without proof. It is, in fact, more difficult to prove than Theorem 2.9 because in deforming 
7 to a projection we have to keep Tr[7] = N. The proof was first given in [LEI] (see [BV] for a 
simple proof). 

2.12 THEOREM (Variational principle): IfV is positive then 

E^^[N) : = inf{f (7) | 7 is admissible with Tr[7] = N} 

= inf{£J(7) \ ^ is an admissible projection with Tr[y] = N]}. (2c.33) 

Moreover, if V is strictly positive then any HF minimizer must be a projection. 

Combining Theorem 2.11 and Theorem 2.12 we see that a HF ground state p for a system with 
strictly positive V is not only particle conserving (a = 0) but, in fact, has fixed particle number, 
i.e., p(A/-2) = p{Mf. 

Theorem 2.12 is useful in many cases for obtaining upper bounds to the Hartree-Fock energy, 
E^^ , and hence an upper bound to the true quantum energy, E^ . It allows one to deal conveniently 
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with what is sometimes called the "orthogonahty problem" . Take any matrix 7 satisfying < 7 < 1 
(as an operator) and Tr7 = N, and compute the one-body energy 

(7) = Ti^h = J2 lij hji (2C.34) 

and the two body energy 

-^2(7) = ! X] Vkl;mn{'ymk7nl - Imllnk) (2c.35) 
k,l,m,n 

Then 

< ^HF < ^ (2c.36) 

The important point here is that the matrix 'ymklnl — Imllnk is not, in general, the two-body 
reduced density matrix of any A^-body density matrix unless 7 comes from a Slater determinant, 
i.e., 7ij = Yl^=i Vf'^'j 1 with (p^,...,(f^ being N orthonormal functions. Nevertheless, (2c. 36) 
continues to be true. 

Contrary to the proof of Theorem 2.9 the proof of Lemma 2.10 is unchanged in the canonical 
case (because the 1-pdm 7^ = (1 — t)j + tj', automatically satisfies Tr[7t] = A'' if Tr[7] = N and 
Tr[j']=N). 

2.13 LEMMA: If V is positive and if p is a HF ground state for the Hamiltonian H under 
the constraint p{M) = N then p is a true ground state for the Hamiltonian Hp under the same 
constraint. 

The mean field Hamiltonian has a simpler form when restricted to fixed particle number. 



'■— 5^ ^ij^l^j + 2 5^ ^kl;mn 
i,j kl;mn 



'ykm'^l'^n + Inl^k'^m Iml^k^n InkA^m 



(2c.37) 



which is simply a (number preserving) independent particle Hamiltonian. 

When V is strictly positive we can prove a much stronger result than Lemma 2.13. Namely, 
not only is a HF ground state, p, for H a true ground state for the mean field operator Hp, but it 
is the unique ground state for Hp satisfying the constraint p{M) = N. 

This uniqueness result is equivalent to the striking statement made in the introduction, that 
no degenerate energy levels of Hp (if they exist) can be only partially filled in the HF state. We 
empasize again how contradictory this is to what is taught in elementary chemistry courses. We 
learn in the theory of chemical binding how two unpaired electrons in separate atoms can pair 
up and create a strongly bound molecule. The concept of unpaired electrons relies of course 
essentially on the independent orbital picture of HF theory. What we prove is that there are no 
unpaired electrons in HF theory. 

The obvious question is of course: But what then has happened to the spin degeneracy? The 
answer is, as already discussed, that there may not be spin degeneracy in HF theory, because even 
the spin symmetry can be broken. Likewise for the angular momentum degeneracy; there may not 
be spherical symmetry. 

The following is the theorem that there are no unfilled shells. It was proved by M. Loss and 
the authors [BLL]. Since the proof is very short we repeat it here. 
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2.14 THEOREM: (Shells are always closed.) Assume that the two-body potential V in 
the Hamiltonian H is strictly positive (as in Coulomb case). If p is a HF ground state (a Slater 
determinant state) for H , i.e., a HF minimizer subject to the constraint P{M) = N , then p is the 
unique ground state of the mean field operator Hp satisfying the constraint P{N) = N . 

Proof: From Theorem 2.12 we know that the 1-pdm 7 of p is an A''-dimensionaI projection, 
which we can write 

N 

i=i 

where gi,g2,. . . are N orthonormal eigenvectors of 7. Since p is a ground state for Hp we can 
assume that gi,g2, ■ ■ ■ are eigenvectors of Hp. 

If there is another ground state with the same number of electrons there must be a degenerate 
level containing a vector from 7, say Qn, and at the same time a normalized eigenvector g' which 
is not in 7, i.e., with jg' = 0. We can then define a new A/'-dimensional projection by 

N-l 

l' = W){9'\ + J2\9j)i9j\- (2C.39) 
i=i 

The HF state p' with 1-pdm 7' is then also a ground state for Hp, i.e., p{Hp) = p'{Hp). 
Since 7 is a minimizer for the HF functional we have 



OKSij')- £(7) -iP'iHp)-p{Hp)) 



kl;mn 



(7 - l')mk{l - i)nl - (7 - i)nk{l " l') 



ml 



(2c.40) 



The operator 7 — 7' = |s'Ar)(5Jv| — \g'){g'\- Hence, since V is strictly positive. 



kl:mn 



(7 - 7')mfc(7 - l')nl - (7 - 7')nfc(7 " l'] 



ml 



={9n ^ gN\V\gN <8) 5iv) + {g' ^ g'\V\g' (g) g') - {gn ® g'\V\gN ® g') 

- {g' ® gN\V\g' (g) 5iv) - {gN ® gN\V\gN ^ gN) - {g' ^ g'\V\g' ^ g') 
+ {gN ® g'\V\g' ® gn) + {g' ® gN\V\gN ® g') 
={gN A g'\V\g' A gN) < 0, (2c.41) 

which contradicts (2c.40).B 

To illustrate this result let us consider the simple, but not altogether trivial example of N = 1 
for the atomic Hamiltonian. In this case the HF ground state energy is in fact the true ground 
state energy. It is simply the ground state energy of hydrogen (with nuclear charge Z), i.e, —jZ'^ 
in our units. Owing to the spin degeneracy the ground state is doubly degenerate. We consider 
the spin up ground state and write it -0 = ]). Here we have used the hydrogen ground state 

wave function (p(x) = CZ^/'^ exp{—Z\x\). It is important to realize that the mean field operator 
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corresponding to is not just the hydrogen operator, but rather the operator that acts on a 
general tp' = \ip'{x)a), with a =t or a =J, according to 

H^\if'{x)a) = (^-A j f{yf\x - y\~^dy^ \f'{x)a) 

- S^^ J ipiy)ip'{y)\x - y\-'dy\^'ix)a). (2c.42) 

Notice that is spin dependent — it does not commute with spin rotations. 

The mean field operator and the hydrogen operator agree on their common ground state ip, 

i.e., 



Z 



\x\ 



H^^l:=[-A- — \i^ = -\Z'i^. (2C.43) 



On the other ground state for hydrogen, namely, |(p |) we find, however, 

{if i \H^\ip i) ={ip i I ^-A - 1^ I) + j j ^{xfip{yf\x - y\-^dxdy 

= -\Z'^ + cZ. (2C.44) 

Thus the mean field operator induces a gap in the energy between the two degenerate ground 
states of hydrogen. 

Theorem 2.14 states that this is not only true for the state \ip |), but in fact, the ground state 
for is unique. More importantly, by Theorem 2.14 this is not particular to the case N = 1. 



III. THE GENERALIZED HF THEORY FOR THE HUBBARD MODEL WITH 
ATTRACTIVE INTERACTION 

Ill.a DEFINITIONS 

In this section the generalized HF theory will be applied to the Hubbard model with attractive 
interaction. Our main result (Theorem 3.12) in Sect. IILe) will be that the HF ground state and 
positive temperature Gibbs state are unique, modulo global gauge transformations. First, we recall 
the definition of the Hubbard model. Let A be a finite lattice, i.e., a finite collection of points, and 
let |A| be the number of these points. The one-particle Hilbert space 7i is the 2|A|-dimensional 
space of spinor-valued functions on A, i.e., the set of complex-valued functions on A x {—1,1}. 
The value of such a function at (y, r) (with y G A, r G {—1, 1}) is f{y,r) and the inner product 
is {f\g) = Yly T fiy^ 'T)9{y^ A canonical orthonormal basis in H is given by the delta functions, 
which we denote by |x,(t). Thus, |x,cj) G 7Y is the function fx,a{y,T) = ^y,xST,a- We define our 
complex conjugation in this basis. In the Fock space corresponding to Ti. wc refer to cJ{\x,a)) 
as cj. ^. We often use the abbreviation j:= -|-1 and |:= —1. An orthonormal basis in JT is given 

by {c|:i,CTi ■ ■ ■ c|;jv,crjvl^) I i^i^^i) 7^ {'^31^3) if ^ 7^ j}) which implies that T is 4l^l-dimensional. 
The Hubbard Hamiltonian for negative coupling is 

x,yEX x& 
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The second term in H- is an attractive interaction among the electrons with position de- 
pendent couphng —Ux < 0. (Our notation here and elsewhere is that Ux > 0.) Other authors 
frequently replace the last two factors in (3a.l) by (c^i-c^-f )(c![.|C^|), but we prefer our form (also 
used in [LLM]) because it preserves hole-particle symmetry; if Ux is independent of x the distinc- 
tion is unimportant. In the first term, the |A| x |A| self-adjoint matrix t = {txy}x,yeA is called the 
hopping matrix. Since we shall later include a chemical potential we can and will assume without 
loss of generality that Tr[t] = 0. We do not assume txx = 0, but when t is bipartite (defined below) 
the condition txx = is automatic. We emphasize that up to this point A was an arbitrary set and 
need not have any topological structure. It is the structure of t, linking different points in A, that 
makes the embedding of A into R'' sometimes useful. Indeed, in the original model introduced 
by Hubbard, Kanamori and Gutzwillcr [HJ], [KJ], [GMC], A is a finite cube of lattice points in 
Z'^ and txy := r when x and y are nearest neighbors, and txy = otherwise. The imposition of 
periodic boundary conditions makes A into a d-dimensional discrete torus on which t is translation 
invariant. Let us remark that, except in Sect. Ill.g, we shall neither assume that t is translation in- 
variant nor that Ux is constant. However, connectedness, reality and bipartiteness of t will play an 
important role in our analysis. These notions can be conveniently characterized by means of path. 
A path is an ordered sequence {xi,X2, ■ ■ ■ ,Xn} of points in A such that txiX2)^x2X3) • • • ,ia;„_ix„ 
are all non-vanishing. We will always assume t to be connected, i.e., any two points x,y & A 
are linked by a path {x,xi, . . . ,Xn,y}- We say t is real if t^y G R for all x,y e A. For real t 
self-adjointness implies that txy = tyx = tyx- For every closed path {xi,X2, . . . , x„, xi} the product 
txiX2tx2X3 ■ ■ ■ tx„xi obviously yields a real number provided t is real. Conversely, if tx^x2 ■ ■ ■ tx„xi is 
real for every closed path then t is unitarily equivalent to some real hopping matrix t' = WtW^ 
where W is a gauge transformation (see [LL] Lemma 2.1). Hence, unless there is at least one closed 
path for which txix2 • • • is not real, we may as well assume that t is real. 

The matrix t is said to be bipartite if there are two disjoint subsets A,BCA with AuB = A, 
such that txy = whenever both x,y £ A or x,y G B. Evidently, t is bipartite if and only if all 
paths contain an even number of points. 

One important property of a bipartite t is that it is unitarily equivalent to —t. Indeed, defining 
the unitary |A| x |A|-matrix (-1)^ = [(-l)'']"^ = [(-l)""]^ by 

(-l)-:=5;|a;)(a;|-5^|a;)(a;|, (3a.2) 
xeA xeB 

one easily checks that 

(-l)^t(-l)^ = -t. (3a.3) 

If it is unambigously clear from the context what is meant, we will use the same symbol ( — 1)^ to 
denote the function which takes the value 1 on the A sublattice and —1 on the B sublattice. The 
Dirac notation is used in (3a.2), whereby |¥')(<p| denotes the projection onto a normalized vector 
ip. The vector x is here the delta function 5x,- at x G A in the space of complex-valued functions 
Ha on A (not A x {T, j}). Let us now make a few remarks about the different Hilbert spaces we 
will encounter in what follows. We start with Ha. Mathematically we may view the one-particle 
Hilbert space H as 

H = Ha e = (g) >iA (3a.4) 
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and, consequently, the space on which the 1-pdm are defined as 

H®H = C^ ^Ha- (3a.5) 

(We wrote equahties in (3a. 4) and (3a. 5) despite of our awareness that an isomorphism would have 
been mathematically more appropriate.) It is thus clear that operators on l~i — ® can be 
written as 2 x 2-matrices with operators on Ha as entries and, likewise, operators on 7^ © H as 
4 X 4-matrices with operators on Ha as entries. The usual conventions for matrix algebras are 
understood and, in particular. 



/ mil 


mi2 \ 




\m2i 


m22 ) 





qmi2 



(3a.6) 



where q and niij are complex numbers or operators on Ha, and likewise for 4 x 4-matrices. It 
will often be convenient for us to change between the Hilbert spaces Ha, H, and H ®H, and to 
take traces over these various spaces. To simplify notations we shall use a common symbol, Tr, for 
these traces; the Hilbert space in question will be evident from the operator whose trace is being 
computed. Likewise, we shall denote the identity operator on these different spaces be the common 
symbol 1. It will always be clear from the context which identity we are refering to. 

We are now in a position to write down the pressure functional for the Hubbard model. We 
choose to define the 1-pdm in terms of the orthonormal basis of delta functions \x,a). Hence, the 

^ j _ 1 < 1 is given by 

S{T) = TVN - J2 ^Aii^ T b\x T) - i |7|x i) - i] 

xeA (3a.7) 

-|(x,T|7k,i)|' + KxTHx 

Here, we identified t with the operator ( ^ ^] onH. 



^0 t ^ 

The thermodynamic pressure of F (actually, the pressure multiplied by the "volume") is given 

by (with P = 1/fcBoltzmannr) 

-VpA^) = s{r) - (3-^s{r) - f,N{r) 

= S{r) + |/?-^Tr[r In r + (1 - r) ln(l - r)] - nTr[j]. (3a.8) 

We have here introduced the chemical potential We could equivalently have replaced the Hamil- 
tonian iT_ by — ^J\f, where J\f is the particle number operator. Also we have introduced the 
notation N{T) := Tr[7] for the particle number expectation in the state described by F. 

Our aim is to characterize the set of maximizing l-pdm's for the pressure (which we typically 
denote by Fq) as explicitly as possible and to determine the pressure V{I3,^) of the system, i.e.. 



P(/3,//) := max{P^,^(F) | < F = 1 ! < 1) = 'P^A'^^ 



(3a.9) 



A priori, the max in (3a.9) should be replaced by supremurnhvA,, since the underlying Hilbert space 
is finite dimensional, the existence of at least one maximizing 1-pdm is assured. 
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Ill.b LINEARIZATION OF THE PRESSURE FUNCTIONAL 

As a first illustration of the notation introduced above, a 1-pdm T as an operator on C^'SiHa, 
is written as 

/7t 7* \ 



1- 



t 



71 



aj -a* 1 - 7t 7* 



it 



(36.1) 



where {x\-f^\y) := {x , a\-f\y , a) , (xlj^ly) := (x, T |7|y,i)! (a^laalz/) := {x,a\a\y,a), and 

(x|Q;*|y) := {x, t |a|y, i). Now, observe that S{T) depends neither on nor on Q!|, and that ^^(r) 

would be lowered, if we were allowed to replace 7* by 0. Indeed, as the following lemma shows, 

a restriction of our attention to the 1-pdm with ai,ai and 7* all equal to zero and, moreover, 

7l = 7| and a* = is justified. Such matrices are of the form (with empty spaces denoting 

zeros) 



a' \ 



7 



-a 



1-7' 



Va't 



(36.2) 



1-7'/ 



where 7' = 7'^ and a''^ = a' . 

3.1 LEMMA: For all /3 > (including the /? = 00) and all /j, we have 

-p(/5, ^) = min {-J'/3,;^(r) | F of the form (3b.2) and satisfying < F < 1 } . 

Proof: Define two orthogonal projections on H ® 7^ by 



P := 



0/ 



(36.3) 



Clearly, PP = and P + P = 1. Observe that for any 1-pdm F written as in (3b. 1), we have 



/7T 



PTP + PVP 



7i 
— al 



t 



—a* 
1 -7t 



(36.4) 



This operator F is also of the desired form (3b. 1). Moreover, < F = PTP + PTP < P + P = X 
and hence F is a 1-pdm. As remarked above, £-{T) < £(T) and A^(F) = A^(F). It remains to 
show that — ^(r) < —S{T). We recall that A Tr{/(yl)} is a concave function of the self-adjoint 
operator ^4 if / : R — > R is concave, i.e., /(Ax + (1 — X)y) > A/(x) + (1 — A)/(y) for all x,y and 



29 



VBEHLJPS - 25/N0V/93 



all < A < 1. Applying this to f{x) := x\nx + {1 — x) ln(l — x), we observe that -S'(r) is concave 
in T; 

5(|ri + iPs) > i5(ri) + i5(r2), (36.5) 

for all 1-pdm Ti and Ts. We evaluate (3b.5) on Ti := T and T2 := (P - P)T{P - P), the latter 
being a 1-pdm since P — P is unitary and preserves the form (3b. 1). But, using the unitarity of 
P — P again, (3b. 5) yields 



s{r) = 5(iri + ir^) > is(r) + i5((p - p)t{p - p)) = s{t). 

Hence, we may restrict the variation to l-pdm T of the form (3b.4). 
Moreover, defining the unitary 



(36.6) 



W 



-1 

V 



1 
-1 0/ 



(36.7) 



we observe that H^FH^^ coincides with F except that and are interchanged as are a* and . 
Hence, 



r := ir + ^WTW^ 



7 
-a't 



—a 

1 -y 



(36.8) 



1-7'/ 



where 7' := |(7| + 7^) and a' := |(q;* + a^) 
As before, concavity of S implies 



s{f) = s(ir + lwTw'<) > is(r) + ls{wrw^) = s{t), 

with equality if and only if f = T. We have Ar(r) = N{f) and 

^(r) - £{f) = U4l[{x\j^\x) + (x|7jx) - If 

X 



(36.9) 



(36.10) 



with equality if and only if (a;|7-|'|a;) = (x|7||a;) for all x G A. ■ 

Thanks to Lemma 3.1 we may now restrict ourselves to 1-pdm of the form (3b. 2) for which 
the 4 X 4-matrix formalism is clearly redundant. Indeed, introducing the unitary operator 



Y :-- 



^10 0\ 

1 

10 

\o -1 o/ 



(36.11) 



otlTKBTC, one easily checks that F in (3b. 2) becomes 



7' 



a 



1-7' 



7 



a' 
1-7'. 



r' 



r' 



(36.12) 
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Although it seems that the conditions that F be a 1-pdm onTi^Tl arc equivalent to the conditions 
for r' to be a 1-pdni on Ha © Ha, there is one important difference. We must require V to satisfy 
< r' < 1 and to be of the form 

with 

(y)t = y , ^ (3^_;L4) 

Because of this difference we will distinguish between Ha © Ha and even though they are 
isomorphic. Physically, H is the space of spin-up and spin-down particles whereas the action of Y 
shows that Ha® Ha, on which V is defined, is rather the space of spin-up particles and spin-down 
holes. 

We shall denote the projection operator in Ha © Ha which projects onto (functions non- 
vanishing only at) the site a; G A by la;. The diagonal part of a 1-pdm V is denoted by 

/y(x) a'ix) \ 
T' := l^r'l^ = 1^, (36.15) 

\a'{x) l-i{x)) 

where ^'{x) := (a;|7|a;) and a'{x) := {x\a'\x). Note that the product of any two operators on 
Ha® Ha of the form A{x)lj; and B{x)lj;, where A and B are 2x2 matrices as in (3b. 15), is given 
by the operator {A{x)B{x)) 1^. According to this rule one easily checks that the trace on Ha®'Ha 
of satisfies 



X 



iTr[r^2]=y(x)2 + |a'(x)p + i-y(x) 



= [y(x)-ip + K(x)p + i. ^^^-'^^ 

By means of (3b. 16) we may rewrite the pressure functional (3a.8) as 

-r^A^') := -VpA^) =2TVN -hJ2Ux {TrK'] + 1) 

X 

- 2/xTr[y] + /?-iTr[r'lnr' + (1 - T') ln(l - F')]. (36.17) 

Finally, we embed the hopping matrix into the 2 x 2-matrix formalism. We define for any real 
number A 

With this definition we obtain 

lTr[TxT'] = TV[iy] - AIV[y] + iA|A|, (36.19) 

and, hence, 

-VpA^') =Tr[T,T'] - i 5^C7.TV[r;2] 

X 

+ p-^Tr[T' In T' + (1 - F') ln(l - F')] + i ^ - ;^|A|. (36.20) 
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We rewrite T' as 



r. = (7'(x)-|) _i Ix+U'r.^ "^^Mlx + |l.. (36.21) 



The cross terms in T'^^ are traceless, which impUes that 

Tr[r;2]=Tr[(r;-iuf] + iTV[U] 

= TV[(r;-il.)2] + l. (36.22) 

Using (3b.22), (3b.20) becomes 

-rf,,,{r') = TT[T,r'] - \ u.mK - 

X 

- n\k\ + /?-iTr[r' Inr' + (1 - r') ln(l - r')]. (36.23) 
We associate a multipUcation operator UD := Y^UxD^lj; with any real function d{x) and 

X 

complex function S{x) by 

fd{x) 5{x) \ 

D^=\ . (36.24) 

\six) -d{x)) 

Notice that has the same form as — ^Ij;. The operator UD is really a matrix- valued potential 
which will enable us to linearize the quadratic trace in (3b. 23) by means of the identity 

-Tr[(r^ - il.)2] = min { - 2Tv[D^{r', - il.)] + Tt[DI]} 

a,o 

= min { - 2Tr[D^T'J + Tr[Dl]}. (36.25) 

a,o 

Indeed it is possible to convert the variation over 1-pdm's in (3a.9) into a variation over all matrix- 
valued potentials D, which we will show by means of the following lemma. 

3.2 LEMMA: Let D be any matrix-valued potential as in (3b. 24), Q be a self-adjoint operator 
on Ha ffi Ha and V the unitary given by 



^-1 0^ 

Suppose that F is an odd real-valued function, i.e., F[x] = —F[—x]. Then F defined by 

f-^:=F[Ti,-UD] + Q- VQV'' , (36.27) 

fulfills (3b. 13), (3b. 14). 

Proof: It is easily checked that FT^Ft = and VUDV'> = -UD. Then by the spectral 
theorem, 



V{r - |)Ft =F [V{T^ - UD)V^ + VQV'' - Q 
=F [-(T^-UD)] - {Q - VQVi) 



- {F [T^ -UD]+Q- VQV^) 

- (r^), (36.28) 
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which is equivalent to (3b. 13) and (3b.l4). ■ 

3.3 THEOREM (Positive temperature pressure): For all < (3 < 00 and all /j, we can 

write the pressure V{P,IJ-) as the following variation over the functions d and 6. 

-V{(3,fi) = min7^;3,^(d, S) - {20-^ \n2 + n) |A|, (36.29) 

d,5 

where 

np,^id,6) :=nfs,f,{D) :=-p-^Tr[lncosh§{T^-UD)] + ^J2u^Tr[Dl]. (36.30) 



7/ a potential D minimizes TZis^n then the operator 



r' = 1 + exp 



P{T^ - UD) 



(36.31) 



is a minimizer for (3b. 23), (i.e., defines a HF Gibbs state) and satisfies the consistency equation 

r; = {D, + i) 1,. (36.32) 

Conversely, ifV is a minimizer for (3b. 23) then the potential D defined by (3b. 32) minimizes Tip ^^j, 
and satisfies (3b. 31). 

Proof: By means of (3b. 25) we write 



-P(/3,/x) = minmin|Tr[(r^ - UD)r'] + i i7,Tr[i^2] 

- At| A| + /J-^Tr [r' In T' + (1 - F') ln(l - T')] | 

= mini min|Tr[(T^ - UD)T'] + /3-^Tr[r' InT' + (1 - r')ln(l - r')]| 
d,d I r' L J 

+ - H\A\. (36.33) 

X ) 



The minimum is over all < F' < 1 satisfying (3b. 13) and (3b. 14). Note that all we did in (3b. 33) 
was to interchange the two minimizations this is of course allowed since we are simply looking for 
the minimum in the set of all F', d and 6. 

First relaxing the minimization to be over all < F' < 1, we can explicitly compute the 
minimum over F' in the second line in (3b. 33). We see that the minimum is uniquely achieved 
for the F' defined in (3b. 31). Wc observe, however, that F' — ^ = — |tanh/3(T^ — UD) is of the 
form suitable for Lemma 3.2 (choosing F[x] = — ^tanh[/3a;] and Q = 0). The operator F' defined 
by (3b. 31) therefore automatically satisfies (3b. 13) and (3b. 14) and we have, indeed, found the 
minimizer in the second line in (3b. 33). Moreover, a simple computation then shows that the right 
side of (3b.33) is identical to the left side of (3b. 29). 
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If D minimizes T^js^^i and we define V by (3b. 31) we have from (3b. 25) and (3b. 29) that 
-V{15, n) = n^AD) - In 2 + //)|A| =Tr[(r^ - UD)T'\ + \ f^xTr[Z)^] 

X 

- n\K\ + /3-iTr[r' In T' + (1 - V) ln(l - V')] 
>-Vp,^{V'). (36.34) 

Since —V{P,^j.) < — P/3.^(r') we must have equality in (3b. 34) but it follows from (3b. 25) that this 
can only happen if (3b. 32) holds. 

Conversely, if F' minimizes —T'/3,n and we define D by (3b. 32) we see that 

X 

- //| A| + ^-^Tr[r' In r' + (1 - r') ln(l - r')] 
>7e;3,^(D)-(2/3-Mn2 + //)|A| (36.35) 

and we conclude that D is a minimizer for TZp^^ and that (3b. 31) must hold. ■ 

We remark that the minimum of Tlis^^{d,6) will be attained for d and S satisfying 

d{x)^ + \6{x)\'^ < I for all a; G A. (36.36) 

In fact, since < F' < 1 we see from (3b. 32) that -D^ — i which implies (3b. 36). 

III.c GAP AND ZERO TEMPERATURE LIMIT 

Our discussion will turn in this section to the zero temperature limit /3 ^ 00. Recall that the 
generalized HF energy is given by 

E^^{n) = mf{e{r')- nN{r') \ r'is a l-pdm}. (3c.l) 

A minimizing l-pdm Fq, that is one that satisfies, -E"^(/x) = £ {F'q) - /iN {T'^) , is called a HF ground 
state. To make contact with our previous notation let us denote 

-V^A^') := ^(F') - fxN{T') = hm -V(,A^'), (3c.2) 

p—>oo 

and -V{oo,fi) := E^^{n) = lim^^^o -V{P,n). 

We can derive an analog of Theorem 3.3 for the zero temperature pressure by simply dropping 
the term — /3~^S'(r), and essentially repeating the whole positive temperature discussion. There is, 
however, one subtle difference. Given the potential D, the minimizing F' in (3b. 33) was uniquely 
defined by (3b. 31). In the zero temperature case the F' we are looking for is simply the minimizer 
of Tr[(T^ — UD)T']. If the operator (T^ — UD) has zero eigenvalues this minimizer would not be 
uniquely defined. In our case, however, we prove in Lemma 3.5 below that (T^ — UD) has no zero 
eigenvalues. In fact, 

|e,| > imin{C/^}|A|-i > 

X 

for any eigenvalue ej of — UD. We can therefore write that the minimizer of Tr[(r^ — UD)r'] 
is F' = x{Tu. - UD) = lim^^oo (1 + exp[/?(r^ - UD)])~^, where x(a) := 1 if a < and zero 
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otherwise. The point is that the value is unimportant. Let us note that there is another way 
of characterizing the zero temperature states, thereby avoiding the repetition of the whole positive 
temperature discussion. Namely, the main theorem 3.12 of this chapter allows us to obtain all zero 
temperature states as limits of positive temperature states. 

3.4 THEOREM (Ground state pressure): We have 

-7^(00,//) = eH'^I/x) = min7^oo,^(c^,^) - /x|A|, 

d,5 

where 

X 

If a potential D minimizes Tlc>o,ii then the operator (with x defined as above) 

T' = x{T^ - UD), (3c.4) 

minimizes (3c.2), (i.e., defines a HF ground state) and satisfies the consistency equation 

K = {D^ + \)lx. (3C.5) 

Conversely, if V minimizes (3c.2) then the potential D defined by (Sc. 5) minimizes Tloo,ti and 
satisfies (3c. 4) ■ 

The reader may wonder why the pressure depends on the absolute value oiT^ — UD and not 
only on the negative eigenvalues. In this context it should be kept in mind that because of the 
special form of the operator — UD the trace — ^Tr|r^ — UD\ is, in fact, equal to the sum of the 
negative eigenvalues. In comparing (3c. 3) to (3b. 30) wc notice that lim^^oo Z?""'^ In cosh(/?x/2) = 
^\x\ it is therefore natural to write the absolute value in (3c. 3). It remains to prove the absence of 
zero eigenvalues. 

3.5 LEMMA (Gap estimate): Let D be a minimizing matrix-valued potential for TZj3^^ 
with < /3 < 00 and denote the eigenvalues of — UD by ei, 62, • . . , e2|A| ■ Then, for any 
j = 1, 2, . . . , 2|A|, we have 

|ej| > iC/„.i„|A|-^ -/?-^21n2, (3c.6) 

where Umin '■= min{i7a;}. 

Proof: If D minimizes '7^/3,^ for finite (3 we know that the operator V defined in (3b. 31) 
minimizes the functional —V/s^^ in (3b. 23). Likewise, if D minimizes TZoo.^ we know that the 
operator F' = lim/3^oo(l + exp[/3(rj^ — UD)])~^ minimizes the energy functional (3c. 2). Notice 
that these definitions of V imply (3b. 13) and (3b. 14). What we do not know a-priori is that if 
P = 00 (i.e., zero temperature), then F' agrees with x{T^ ~ UD). 

Let Q := \{p){{p\ denote the projection onto some eigenvector ip oi — UD with correspond- 
ing eigenvalue e. Note that VQV^ is the projection onto the eigenvector Vlp of — UD with 
eigenvalue — e, provided V is the unitary given by (3b. 26). 

In terms of the operators F' discussed above we define 

f ' := T' + 6Q- dVQV^ . (3c.7) 



35 



VBEHLJPS - 25/N0V/93 



We choose -1/2 < S < 1/2 such that eS > 0. With this choice < T' < 1. Note that (p is 
an eigenvector of T' . We denote its eigenvalue by A. Then (p is also an eigenvector of T' with 
eigenvalue A + c^. With our choice of 6 we have 0<A, A + 5<1. 

Lemma 3.2 now implies that T" is admissible on Ha © Ha, as it corresponds to the 1-pdm 



Since F' minimizes —Vjs,^ it follows that 

= - 5Tt[T^{Q - VQV^)] - \ ?7,Tr[(r; - ^l^f - (r^ - il, + 6Q^ - SVQ,V'<f] 

X 

+ 2/?-^[Aln A + (1 - A) ln(l - A) - (A + (5) ln(A + (5) - (1 - A - 5) ln(l - A - 5)] 
= -2Se+ i J2 Yl UxTr[{Q, - VQ.V^f] + 20-' [/(A) - /(A + S)], (3c.9) 

X 

where /(A) := A In A + (1 - A)ln(l - A). Clearly, -ln2 < /(A) < for < A < 1. Hence 
|/(A) — /(A + (5)1 < ln2. We insert this estimate into (3c. 9) and obtain, using 5e = \d\\e\ 

2\e\ > l\6\Y,UxTi[iQ,^-VQ,V''f]-2\5\-'f3-'lu2. (3c.l0) 

X 

Furthermore, 

IV[(g, - VQ^V^f] = 2 T)P + Mx, (3c.ll) 

leads us by the Schwarz inequality to 

YUxTriiQx - VQ^V^f] >2C/^i„5^ (|</p(a;,T)|' + I^^C^'i)!')' 

X X 

>2C/„i„|A|-' ( j;iv^(x,T)|' + k(x,i)p 

\ X 

=2U^i^\A\-\ (3C.12) 

Inserting (3c. 12) into (3c. 10) and choosing \S\ = | concludes the proof. ■ 

We interpret any minimizing F in Theorems 3.3 and 3.4 as the physical state of the system at 
(3, n. Since F, subject to (3b. 31), is the Fermi distribution of the corresponding operator — C/D, 
we might as well interpret its eigenstates as the only orbitals the electrons can possibly occupy. This 
is nothing but viewing T^ — UD as the relevant quasiparticle Hamiltonian for the considered system. 
From the BCS theory of superconductivity the question arises whether or not the quasiparticle 
spectrum, i. e., the spectrum of — UD has a gap around for low enough temperatures. The 
answer to this question is always positive as we have just pointed out in Lemma 3.5. We remark, 
however, that our estimate for the size of the gap becomes trivial in the thermodynamic limit. In 
Theorem 3.15 below we give a formula for the gap in the translation invariant case. 

Ill.d BROKEN GAUGE SYMMETRIES 
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We pause to remark that so far no assumption was made about the hopping matrix t, except 
that it be self-adjoint, connected and tracclcss. In Theorems 3.3 and 3.4 in the previous section 
we estabhshed a unique correspondence between the HF ground and Gibbs states (described by 
r') and the potentials D that minmize TZp^^. We can therefore now entirely concentrate on the 
determination of the functions d and 6 that yield a minimizer D for 7?.^^^, where < /? < 00. 

A key ingredient in the following analysis are the representations 



In cosh X = In 



fc=0 



1 



vr(fc + 1/2) 



00 / ^2 



In 1 + ^ Uc 



(3c/. 1) 



The significance of (3d.l), as examples of integrated Pick functions, was emphasized in Lieb and 
Loss [LL] and Kennedy and Lieb [KL]. The virtue of (3d.l) is, roughly speaking, that it allows us 
to convert traces into determinants. In fact, using these representations we may write 

00 

n^AD) = - 5^[lnDet{47r^(A; + l/2fr^ + (T^ - UDf} 

k=0 

- 2|A| ln{47r2(fc + 1/2)2/3"'}] + ^ E U.Tr[Dl], {3d.2) 

X 

n^,^{D)=-^J^ [i In Det{c' + (r^-C/D)'}-2|A| In c](ic 

+ i5]C/,Tr[I)2]. (3d.3) 

X 

Since we are dealing with finite dimensional matrices, convergence of the above expressions are 
evident for any choice of D. 

3.6 LEMMA (Phase alignment): Lett be real. Then, 

nfsAd,S)>nfs,^{d,e''\6\). (SdA) 

If S{x) 7^ for all x € A, equality holds in (3d.4) only if S{x) = |(5(a;)|e*^ for all points x e A and 
for some fixed < 9 < 2tt. 

Remark: Lemma 3.8 below will show that, for the ground state, either 5{x) ^ for all x or 
else S{x) = for all x. 

Proof: In order to prove (3d.4) we may assume that d{x) ^ for all x G A. In fact, if we have 
proved (3d.4) for non-vanishing 6 we conclude for general S by the continuity of TZp^/j, that 

Hf^Ad^S) = lim7^;3,^(c^,5 + e) > lim 7^;3,^(ci, |5 + e|) = 7^^,^(ci, \S\). (Sd.b) 

For the cases of equality, however, we have to assume 5 non- vanishing. 

Just as D corresponds to d and S via (3b. 24), we denote the operator corresponding to d = d 
and 5 := e*^|(5| by D. In view of (3d. 2-3), it suffices to prove, for every real c, that 

InDetfc^ + (T^ - UDf] < InDetfc^ + (T^ - UDf] (3ci.6) 
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and to show that equahty in (3d. 6) implies S{x) = e*^|(5(a;)| for some < ^ < 27r. 
To this end, we define an operator 

X 

with Bx > (to be chosen below). Furthermore, we rewrite UD in the following way. Let 
F := Ex -^xlx and G := J2x ^x^x, where 



Fx := Ux d{x 
Gx '■= Ux 

Hence, 



1 

d{x) 



T^^UD = Tf, - F - G. (3(1.9) 

Notice that both F and G commute with B, and that F and G anticommute. Thus, abbreviating 
B^/'^QB^^'^ =: Q for any operator Q we have the following identity. 

lnDet[c^ + (T^ - UDf] + 21nDet[B] 

= hiJyei[B^'^{ic + T^-F- G)B{-ic + T^-F- G)B^/^] (3d.l0) 
= \nBet[A{c,D)-{f^,G}], 

where {A, B} := AB + BA is the anticommutator and where 

A{c,D) := c'B^ + (f^ - Ff + G^ + ic[B,f^]. (3d.ll) 

We shall now prove that 

lnDet[c2 + (T^ - UD)^] < lnDet[^(c, £>)] - 21nDet[S], (3d.l2) 

and that equality holds in (3d. 12) if and only if {T^, G} vanishes. 

We will prove this claim by using a concavity argument. Introducing a unitary 

V--={1 ~o)^ (3d.l3) 

one easily checks that VBV^ = B, Vf^Vt = -f^ = -f^, using the reality of t, and VFV^ = -F, 
YGV"^ = G. Hence, we obtain 



VA{c, D)V^ = A{-c, D) = A{-c, D), V{T^, G}V^ = -{T^, G} (3d.l4) 

Now, from (3d. 10) we see that the determinant in question is real and depends only on c^. There- 
fore, 

lnDet[A(c, D) - {f^, G}] = lnDet[^(-c, D) - {f^, G}] 

= lnDet[^(-c,L>)-{f^,G}] = lnBet[V {A{c, D) + {f^,G})V^] 

= \nDet[A{c,D) + {f^,G}]. (3d.l5) 
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By the strict concavity of lnDet[-], i.e., |lnDet[a] + ^lnDet[6] < lnDet[^(a + b)], with strict 
inequahty unless a = b, we arrive at the assertion in (3d. 12). 

We now choose = U~^\6{x)\~^ , keeping in mind our assumptions that 6{x) 7^ and > 0. 
The key observation is that 

is unitary for all x G A, i.e., B^G^ = 1. Thus, for any two points x,y G A the 2x2 matrices 
belonging to (G)^ and {G)y are identical modulo a phase factor S{x)S{y)\d{x)S{y)\~^ . The last 
term {T^, G} vanishes if and only if 

t^yS{x)\d{x)\-^ = Uy5{y)\S{y)\-^ {3d.l7) 

for all x,y e A. It is essential for (3d. 17) that t be real. From (3d. 17) we learn by the connectedness 
of t that T^G + GT^ vanishes if and only if 5 is of the desired form, S{x) = e^^\S{x)\. 

Notice now that A{c,D) = A{c,D) and that for D we have equality in (3d. 12). Inequality 
(3d. 6) is therefore a consequence of (3d. 12). ■ 

We remark that by (3d. 12) we also get a lower bound on 7^^^^(D) that depends only on 
d and \S\, even in cases in which t is not real. However, this lower bound cannot then be ex- 
pressed as TZp^^{D') for some matrix-valued potential D' . The reason is that (3d. 17) now reads 
txy^{x)\5{x)\~^ = txy5{y)\8{y)\~^ . This equation cannot be satisfied unless (5 = as shown in 
(3d.26) below. 

The strategy for proving Lemma 3.6 can actually be applied to other types of hopping matrices, 
as the following lemma shows. Using the notion of pseudo-spin discussed below (in (3d. 50-53)) the 
conclusion of the next lemma is that minimizers will have aligned pseudo-spin. 

3.7 LEMMA (Pseudo-spin alignment): Let t be bipartite (but not necessarily real) and 
II = 0. Then 

7^^,o(d, 6) > n(3,o [i-ir^/d^+W, 0) . (3d.l8) 



In case t is not real and n(x) := y^P(x) + |5(x)P 7^ for all x E A equality holds in (3d. 18) if 
and only if 5{x) = and d{x) = (— l)^n(x) or d{x) = (— l)^n(x) for all x G A. In case t is real 
and n{x) 7^ for all x E A equality holds in (3d. 18) if and only if for all x E X 

(d{x) 6{x) \ ({-irn{x) \ 

\=wA 4, 3d.l9 

\6{x) -d{x) V -{-lYn{x) 



{-ly \ ( {-ly 



where 



for some unitary 2 x 2-matrix w independent o/x G A. 

Proof: Since the main idea of the proof is the same as that in Lemma 3.6, we shall use the 
notation therein. In proving (3d. 18) we can assume n{x) non- vanishing for all x, otherwise we use 
a continuity argument as in Lemma 3.6. For simplicity we write Tf^=o = T. Again, our assertion 
follows by showing that for all real numbers c the inequality 

InDetfc^ + (T - UDf] > InDetfc^ + (T - UD)^] (3d.20) 
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holds, and there is equahty if and only if d and S fulfill (3d. 19) and D is assumed to be some 
matrix- valued potential for which the corresponding functions d and 6 do obey (3d. 19). 
Again, we define Q := B^^'^QB^^^. In analogy with (3d. 10) we obtain 

InDetfc^ + (T - UDf] + 21nDet[B] 

= In Bet[B^/^{ic + T- UD)B{-ic + T- UD)B^/^] (3d.21) 
= ]xiDet[A,{c,D)-{f, UD}], 

where now 

A^{c,D) := {UDf + c^B^ + ic[B,f]. {3d.22) 

Thanks to our assumption of bipartiteness of t and the special choice of ^ = for all x € A, we 
can now repeat the concavity argument above, using the unitary (— 1)"^. Indeed, (— l)^r(— 1)^ = 
-f, {-1)^B{-1)^ = B, {-1)^UD{-1)'= = UD, (-l)^y4*(-c,L')(-l)^ = y4*(c,L'), and thus 

lnDet[^*(c,L»)-{f,J7i)}] = lnBet[A^{-c, D) - {f,UD}] = lnBet[A^{c, D) + {f,UD}], (3d.23) 

Hence, concavity of lnDet[-] again implies that 

In Det[c2 + (T - UD)^] < In Det[^* (c, £>)] - 2 In Bet[B] (3d.24) 

with equahty if and only if {T,UD} = 0. We now make the choice B^ = Ux~^n{x)~^. The 
condition {T, U t)} = becomes 

d{x) /n{x) = —d{y) /n{y) for x G A,y ^ B, {3d.25a) 

txyd{x)/n{x) = t^S{y)/n{y) for x G A, y G B. {3d.25b) 

But (3d. 25) is equivalent to (3d. 19) provided t is real. Conversely, in case t is not real, as we pointed 
out in the definition of the hopping matrix, there exists at least one closed path {x,xi, . . . Xn,x} 
such that txxi ■ ■ ■ txnx is not real. Thus, iterating (3d. 25b) along this closed path we obtain 

!"^^---S^^ ^ = ^ , (3d.26) 

ixxi • • • txnXi "fi^x) n{^x) 

whose only solution is d{x) = 0, implying (5 = on the entire lattice A because t is connected. 
Finally, we remark that we may achieve the bound in (3d. 20) by choosing d{x) = (— l)^n(x) and 
6ix) = 0. ■ 

In order to use Lemma 3.6 to conclude that a potential D minimizing TZfi^n is phase aligned 
we must show that for such a potential 6{x) ^ for all x (or otherwise vanishes everywhere). 
Likewise, in order to use Lemma 3.7 we must show that a minimizing potential has n{x) 7^ for 
all X. We prove these two results in parts (a) and (b) of the next lemma, again using the method 
of choosing an appropriate normalizcr, B. 

3.8 LEMMA (Nonvanishing of minimizers): (a) Ift is real and if D = {d,6) minimizes 
T^I3,ti then either 6{x) for all x or 6{x) = for all x. 

(b) If t — ^ is bipartite (but not necessarily real) and if D = {d,S) minimizes T^is,^ then either 
n{x) ^ for all x or n{x) = for all x. 
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Proof: (a) Let D = {d, S) be a potential such that both sets 

Ao = {x G A I 6{x) = 0} and A \ Aq = {x G A | 5ix) ^ 0} 

are non-empty. We shall prove that D cannot be minimizing for T^p,^- For < r < 1 we define a 
new potential D^- by 

dr{x) = d{x) and ,5,(x) = |^/-^'l'^(^)l' , (3d.27) 

I T(!)av, X G Ao 

where = (Z^xsAq ^ ^ t^x|<^(3:^)P > 0. We shall show that for r small enough 1Zi3^^{D) > 
According to (3d. 2-3) we can write 

POO 

T^fsAD) = - / lnDet[c2 + (T^ - UDf] - 4|A| Inc d/x^(c) + V U^Tr[Dl], (3d.28) 

Jo J X 

00 

where the measure d/ifj for /3 < 00 is a sum of Dirac delta functions diJ,j3{c) = l3~^S[{2k + 1) 

fc=o 

Tr/?"-*- — c]dc, while for /3 = 00 it is d/^oo = (27r)~-^dc. 

We chose S^v to make the last term in TZp.^ independent of r, i.e., [/^Tr[L'^] = J2x ^x'^A^t,x]- 
To conclude that TZp^^{D) > 7^^^^(I?^) we must therefore show that for r small enough 

00 

J {lnDet[c2 + (T^ - UDf] - lnDet[c2 + (T^ - C/L>^)2]}d//^(c) < (3d.29) 


To prove (3d. 29) we again appeal to (3d. 12). We choose = U~^[6r{x)\~^ and obtain that 

InDetfc^ + (T^ - UDf] < lnDet[y4(c, D)] - 21nDet[S] (3d.30) 

and 

InDetfc^ + (T^ - UD^f] = lnDet[A(c, L>^)] - 21nDet[S]. (3d.31) 

Since {T^, Gr} = (here G-r is defined as in (3d. 8) but with S replaced by Sr) we see from (3d. 10) 
that 

A{c, Dr) = (icB + - UDr){-ic + T^- UD^) 
and A{c,D) = A{c,Dr) - 1 + Therefore (3d.30) and (3d.31) imply 

lnDet[c2 + (T^ - UDf] - InDetfc^ + (T^ - UD^f] 

< lnDet[l - {icB + f^ - UDr)-\l - G^)i-icB + f^ - UDr)'^] {3d.32) 
= Trln[l - {icB + f ^ - UDr)~^{l - G^){-icB + f ^ - UDr)~^]. 

Using the inequality ln(l — A) < —A we obtain 

lndet[c2 + (T^ - UD)^] - InDetfc^ + (T^ - UDr)^] 

< -Ti[{icB + f^- UDr)-\l - G^){-icB + UDr)-^] 

< -Tr[ii:^(c)], (3d.33) 
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where we have denoted 

Kr{c) :=(ic5 + - UDr)-\l - G^){-icB + UD,)-^ 

=B-^/^{ic + T^ - UDr)~^{B-^ - BG^){-ic + T^ - UDr)~^B-^/^. (3d.34) 

The estimate (3d. 29) follows if we show that for r small enough 

J Tr[Kr{c)]dn0{c) > 0. (3ci.35) 

We denote by Pq the projection onto the sites in Aq, i.e., Pq = Yl and Pq = 1 — Pq = 

xeAo 

^ Ix- If we assume r < | we then have 

x^Ao 

S-1 _BG^=J2 rUJ^A. - -7=^Ux\Kx)\^x > arPo - br^Po, (3d.36) 

where a and b are strictly positive constants depending on the values of \d{x)\ and for all x. 
We therefore find 

TV[i^^=)] >TaTr[B-^/^{ic + T^- UDr)-^Po{-ic + T^- U Dr)'^ B-^'''] 

- THTT[B-^'''{ic + T^- UDr)-^Po{-ic + T^- UDr)-^B-^'\ 

Since a'Po < < b'l for constants a' and b' (again depending only on \S{x)\ and Ux for all x) 
we have the estimate, with = AA^, 

TV[i^|^)] > Taa'Ti[\Po{ic + T^- i7L>0"'^o|'] - THb'Ti[{c' + (T^ - UDrf)-^]. (3d.37) 

Since the problem is finite dimensional it is clear that the eigenvalues of — UD^ converge 

to the eigenvalues of — UD. Hence, the gap estimate (3c. 6) implies that for r small (depending 
on Ux, Dx and txy) there is a constant gp, satisfying 5^ > for /3 large enough (in particular for 
/3 = 00) such that (T^ - UDrf > gj. Therefore 

[c^ + (T^ - UDrf]-' < (c^ + gj)-H. (3d.38) 

Recall now that the measure d/if) appearing in the Pick representation (11.2) is supported 
away from zero when /3 < 00. It therefore follows from (3d. 38) that for small enough r 

j Tr [(c^ + (T^ + UDr)^) d/x^(c) 

is bounded independently of r for all /3. 

In order to conclude (3d. 35) (and hence the lemma) from (3d. 37) it only remains to prove that 

liminf / Ti[\Po{ic + T^ - UDr)~' Po\'^]dii^ic) ^ 0. (3d.39) 

T^O J 
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Since t is connected we know that 

Po{ic + T^- UD)Po = PoToPo + and Po(ic + - VB)Pq = PoTqPo ^ 0. 

Therefore we must have Po{ic + T^j^ — U D)~^ Pq ^ (because if we have a matrix ^ ^ with 
operator valued entries and c 7^ and b ^0 then its inverse ( " ^] must have P ^ and 7 7^ 0) . 



7 S 

Since lim^^o Tr[|Po(«c + T^- UDr)~'^Pof] = Tr[|Po(«c + - C/-D)"^Po|^] we obtain by Fatou's 
Lemma that 

liminf J TT[\Po{ic + T^ - UD^)-' Pol^W^ic) > J Tj:[\Po{ic + T^ - UD)-^Po\^]dn(}{c) > 
and part (a) follows. 

(b) Part (b) is proved in the same way as part (a). We this time choose D^- with Sr{x) = 
and dr{x) = {—l)^nr{x), where 



nr{x) 



s/l — T'^n{x), X ^ Aq 
rnav, a; G Ao 



Here nav is again defined such that Uxn{xY = UxfiT-ixY . The rest of the proof is identical to 
the proof of part (a) except that we use (3d. 24) instead of (3d. 12) and we choose = C/^nT-(a;)~^. ■ 
Wc pause to discuss the symmetry aspects of the Hubbard model and how they are reflected 
by Lemmas 3.6 and 3.7. Independent of the hopping matrix t and the chemical potential /i the 
Hamiltonian H will always be invariant under a global (i.e., all spins are equally transformed) spin 
transformation Wg = Ws(u') 

Ws4,a>Vt:= 5^ (3d.40) 

a' 

for all a; G A, where u; G C/(2) is a unitary 2 x 2-matrix. Note that Wg is a Bogoliubov transfor- 
mation. 

As always we shall be particularly interested in the spin rotations, i.e. transformations Ws{w) 
corresponding to w in SU{2). The full ^7(2) group is generated by the SU{2) subgroup together 

f e'^ \ 

with the subgroup (isomorphic to C/(l)) consisting oi all wq = [ j, < 6 < I-k. The 

Bogoliubov transformation Wgi := Y^s^Wd) is a global phase change, i.e., 

for all a; G A. Spin rotations and phase changes (i.e., the full C/(2) group) exhaust the symmetries 
of the Hamiltonian R unless we assume more about the hopping matrix t. 
The spin rotations are generated by the (quadratic) spin operators 
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The U{1) phase change is generated by the number operator M = a ^x,a^x,a- 

According to (2.a9) the Bogohubov transformation Ws{w) has its counterpart acting on Ti^BTi 
which we wih denote by Ws{w) and in the special case of a phase change by W^. Indeed, 



Ws{w) = ( _ ) and We 




(3d.42) 



e 



i0 j 



From Lemma 3.1 it is easy to conclude that any V minimizing the energy functional (3.37b) 

corresponds to a HF ground state with total spin zero, i.e., a state with p{S\ + 51 + '^3) — 0- 

3.9 THEOREM (Zero total spin): A l-pdm F corresponds to an SU{2) invariant HF 
state if and only if it has the form (3b. 2). In particular, the HF ground states have total spin 
zero. 

Proof: The condition that a state be SU{2) invariant is that F commutes with Ws{w) for 
all w € SU{2). Recalling that the only 2x2 matrices commuting with all elements in SU{2) 
are multiples of the identity matrix it follows easily that exactly the matrices of the form (3b. 2) 
commute with all Ws{w). It then follows that the HF ground states have total spin zero since they 
are pure states. ■ 

The HF states are, however, not necessarily invariant under the U{1) phase symmetries. If we 
use the unitary operator Y defined in (3b. 11) we can write 



YWeY'^ 

where 



W'e 
' ' 



W^=(^'q (3d.43) 

A l-pdm F will therefore commute with Wq if and only if the corresponding F' commutes with 
W'q. In case of the minimizing Fq it follows from (3b. 31) (for positive temperature) or (3c. 4) (for 
the zero temperature ground states) that Fq commutes with W'q if and only if the corresponding 
matrix- valued potential Dq commutes with W'q. Prom (3b. 24) and (3d.43) we find 

(d{x) 5{x) \ , / dix) e'^^Hix)\ 
W^D^W'^ = W'q[ __ ^'\w'^=\ ^'\. (3d.44) 

' \5{x) -d{x)) ' \e--''H{x) -d{x) ) 

Thus Dq commutes with Wq if and only if 5{x) = 0, which not surprisingly is exactly the condition 
that the state is normal. Notice that (3d.44) agrees precisely with the characterization of the 
possible minimizers given in Lemma 3.6. 

Note also that (3d. 44) yields the same l-pdm for 6 + 7r(mod 27r) as for 0, meaning that the 
representation of the group U{1) by the minimizing states is not faithful. This reflects the fact, 
that generalized HF-states p obey the particular restriction p(ei . . . e2k+i) = 0, where is either 
a c or a c^. 

Now, let us consider the case in which t is bipartite and = 0. Whether t is real or not the 
Hamiltonian H is invariant under the particle-hole transformation 

Wpha4,a>Vph := {-lYac^,a, 
>VphQ!C^,aWph := {-ly-acl^^, 
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where a is any complex number which wc insert to indicate that Wph is antiunitary. In fact, for t 
non-real the particle-hole symmetry of the Hamiltonian cannot be unitarily realized. On TC (BTi, 
with the fixed basis {{x,a)\x G A, a =T,i} in H, Wph corresponds to (the 4x4 operator) Wph, 
acting on operators as [A, B, C and D are here 2x2 operators) 

As for the ^7(1) symmetry we would like to write Wph in terms of a symmetry on Ha Ha- 
This is only possible if we first compose it with a spin rotation which after all will leave our states 

invariant by Theorem 3.9. In fact, 'd w = ( ) have 



ph. 

where W^^^ is the antiunitary map on Ha ® Ha with matrix (in the standard basis) given by 

1 



In order to understand the transformation property of a minimizing Fg it is again enough to 
consider the potential. 

/d S \ , ( -d -5\ fd 5 \ 

-J^^-; H-j J^-u -J 

since (-1)^L»(-1)^ = D. 

If t is not real we compare this with our result from Lemma 3.7. According to the condition 
on the minimizing D therein, we must have 5{x) = for all x G A. The minimizing D is, hence, of 
the form D = D^'i-x or D = ^ —D^'^x where 

X X 

Thus, assuming that t is bipartite and non-real and fi = 0, the minimizing l-pdm T has no pairs 
and corresponds to a normal state. It is unaffected by spin rotations and phase changes which 
means it has spin zero and fixed particle number. Only the particle-hole symmetry is broken due 
to the global change d{x) — >■ —d{x) for all x e A. 

In case of a real, bipartite t at /i = the symmetry group of the Hamiltonian is even larger. 
Let us introduce the Bogoliubov transformation Wbp by 

Whr,cKwl =c\, 

f (3d.49) 

^^AiM, = (-i)^c.i- 

The global pseudo-spin rotation Wps is given by 

>Vps(«;) = W^pW,(«;)>Vbp, (3d.50) 
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for any spin rotation Ws{w), with w G SU{2). It leaves the Hamiltonian invariant in the real 
bipartite case at ;U = 0. 

As in the spin case we could of course have considered the full group of U (2) pseudo-spin trans- 



formations. The transformation Wpsiwe) corresponding to we 



the spin rotation corresponding to wq :- 



Q ) is, however, equal to 



, i.e. yVpsiwe) =Ws{we). Conversely, the 



U{1) group of phase change symmetries is really a subgroup of the SU{2) pseudo-spin rotations. 
In fact, the pseudo-spin rotations are generated by the pseudo-spin operators 



a;eA 



Notice that in terms of how they transform operators, i.e., A ^-^ WAW^ we cannot distinguish the 
unitaries 

exp(|A/'6') and exp(i(S3^). 



and is, therefore, contained in the group of 



We also point out that for real t the particle-hole symmetry may be unitarily realized as 

= W3 {w)Wi^W, (^^) with w=(^^^ I 

spin- and pseudo-spin rotations. 

The unitary operator Wps(w) on H (BH corresponding to a w e SU{2) which we write as 



w 



_ _ given by 



tO-|-|- 

\-{-irwn 
Thus,yi^p.Hyt = (^;sH ^^^^^^ 

wM = I 



, where 



J 



(3d.51) 



w 



TT 



(3d.52) 



In other words, Wp^{w)DxWpg{w)^ is the transformation of the potential D^. Now observe that 
if fulfills (3d.l9) then so does W(,^{w)D^W(,^{w)\ indeed, starting from D with 6 = 0, (3d.l9) 
simply states that all minimizers are of the form 



W'iw)DW'iw) 



t 



(3d.53) 



with an arbitrary SU (2)-matrix w. Hence, for real, bipartite t at = the situation is as follows. 
Whenever S is different from 0, the minimizing F does not correspond to a normal state and the 
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particle number is broken. Likewise, the pseudo-spin rotation is a broken symmetry, as can be seen 
from (3d. 53). 

The complete tables to illustrate the broken symmetries can be found in Chapter V. 

IlI.e SPATIAL UNIQUENESS OF MINIMIZERS 

In the previous section we saw that if a minimizing potential D is nonvanishing there is a whole 
family of minimizers related by global gauge transformations of D. In fact, the phase of 5 or the 
pseudo-spin must be aligned over the lattice and we can only allow global gauge transformations. 
The aim in this section is to prove that all degeneracies of the minimizing potentials D are caused 
by these gauge transformations. Since the HF ground and Gibbs states are uniquely determined 
by D it shows that the only degeneracies of these HF states are due to the symmetries discussed 
in the previous section. 

More precisely, we shall show that if Di and D2 are two minimizers for T?-/?,^ then Df = -D|, 
i.e., ni = 77,2, where as before ni,2 = \Jd\^2 + l<^i,2P- We show this in the case where t is bipartite 
and = 0. If t is not bipartite or /x we also get di = 6,2 (see Lemma 3.11). The uniqueness 
statement in the bipartite case with /x = is an immediate consequence of the strict convexity we 
now prove. 

3.10 LEMMA (Convexity of TZ^fi): If t is bipartite (not necessarily real) and /i = then 
^/3,o {[—^Y'Tt'^ 0^ , regarded as a functional of the function rj = is convex. It is strictly convex at 
r] if r]{x) 7^ for all x^ . 

Proof: We shall prove that if no = y/ Anf -I- (1 — A)n2 for some < A < 1 then 



with equality if and only if ni = 71-2 = no- Since the term Tr[D^] is linear in we may ignore 
it here and we only need to consider the first term in (3d. 2-3). Let Dj, j = 0,1,2 denote the 
potentials corresponding to {dj,dj) = ((— l)^nj, 0), for j = 0, 1, 2. 

To prove the strict convexity we assume no (a;) 7^ (i.e., either ni{x) 7^ or n2(a;) 7^ 0) for all 
X G A. If no vanishes somewhere we still get convexity (but not strict) by a continuity argument. 
For both and D2 we use (3d.24) with B = U~^no{x)~^ and we find 



Since D"^ = njl we get from the definition (3d.22) that A^{c,Di) + A^{c,D2) = A^{c,Do). By 
strict concavity of lnDet[ • ] we obtain 



7^;3,o((-l)"no,o) < A7e/3,o((-l)"ni,o) + (1 - A)7^^,o((-l)"n2,o) 



(3e.l) 



InDet + {T-UDjf < lnDet[y4*(c, D^)] - 21nDet[B]. 



(3e.2) 



A InDet + {T-UDif +(1- A) InDet + (T - C/L>2)^ 



< InBet A^{c, Do) - 2 In Det [B] 



InDet c^ + {T-UDof 



(3e.3) 



^ The convexity is, in fact, strict for all r], but this is more complicated to prove. 
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with equality if and only if ni = n2 = uq. The last equality in (3e.3) follows from {T, UlDo} = 0. 
Inserting (3e.3) into (3d.2-d3) gives (3c.l).B 

If we do not assume bipartiteness and = the situation is somewhat more compicated. We 
have the following convexity result. Note that assuming t — iJ, not bipartite includes bipartite t 
with /^ 7^ 0. 

3.11 LEMMA (Convexity of 7^^,^): Ift is realTZp^^id, Vn'^ - cP), regarded as afunctional 
of d and rj = n'^ is convex, but not always strictly convex. If t — jj, is not bipartite the functional 
Tl/3^n{d, \/n^ — d^) is minimized by unique functions d and n. 

Proof: We first prove the convexity. Note that IZjS^fj, as a functional of the functions r] = 
and d is defined on the convex set d) | |dp < r]}. Given dj and rij for j = 1,2 satisfying 



\dj\ < rij. Define do = Xdi + {1 — X)d2 and no = y^An^ + (1 — A)n2. Our aim is to prove that 
('^o, \/nl- dfj < XTZfs^^ (di, ^nf - dj^ + (1 - A)7^;3,^ (d2, 



(3e.4) 



Let Dj for j = 0,1,2 correspond to {dj,dj) = {dj, ^n'j — d'j). Since we can otherwise use a 
continuity argument we may assume that So(x) = ^yno{x)'^ — do{x)^ 7^ for all x E A. 

For both Di and D2 we use (3d.l2) with B = U-^\6o{x)\-^ = U-^ (no{x)^ - do{x)^^ ^ ^ We 
obtain 



InDet 



c^ + {T-UDjf 



< InDet 



A{c,Dj) -21nDet[S], 



(3e.5) 



with equality if and only if {T^, Gj} vanishes. Using A{c,Di) + A(c, -D2) = A{c,Dq) 



A InDet 



+ (T - UDif + (1 - A) InDet + (T - UD2) 



< InDet 



A{c, Do) - 2 In Det[S] = In Det + (T - UDof 



(3e.6) 



The last equality in (3e.6) holds because {T^, Go} vanishes. The convexity in (3e.4) is an immediate 
consequence of (3e.6). 

We have equality in (3e.6) if and only if A{c,Di) = A{c,D2) = A{c,Do), i.e.. 



-{T,,.. Fi} + Gl + = '{T^,F2} + Gl + Fi = -{T^,Fo} + G^ + F^ 



(3e.7) 



To show that 7^/3,^ need not be strictly convex consider the case of bipartite t with ^ 7^ 0. Let 
di{x) = —U~^ii + (— l)^c and d2{x) = — — (— l)^c for some c > and let 5i = 62 = 0. If we 
take A = 1/2 we find do = \di + \d2 = -U~^ii and Uq- dl = \d\ + \d2 - do = c^- 0^ other 
hand by pseudo-spin invariance we see that strict convexity is violated: 

7^;3,^(dl,o) = 7^^,^(d2,o) =7^^,o((-lrc,o) + J^i7-V' 

X 

=7^^,o(o, c) + Y, u-'^i^ = 7^^,^ (do, 

rr. ^ 



nl-dl 



Assume now that D^ and D2 are two minimizers for TZf^^/j, and define Do as above with 
A = 1/2. By convexity of TZf3^f^ wc conclude that Do is also a minimizcr. Moreover, we know from 
Lemma 3.5ab that either 6o{x) = for all x or 6o{x) ^ for all x. If 6o{x) = for all x we have 

= nl-dl = - d2f + ^Sj + ^dl 
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Hence di = 6,2 and 5i = 82 = ^ and thus ni = n2- 

If 8o{x) ^ for all x, we know that (3e.7) is satisfied and that 

{f^,G,} = (3e.8) 

for j = 0, 1, 2. From the off-diagonal part of (3e.7) we conclude that for all a; 7^ y in A with t^y 
we have 

^ix + Ay = Fox + Foy. (3e.9) 

Hence, Ax - ^Ox = -{Fiy - Foy). If the off-diagonal part t' of t (i.e., t'^y = t^y - t^^S^y) is 
not bipartite we can find a path z-i,Z2,-- - , Z2k = z\ with an odd number of points such that 
^zjZi+i 7^ 0- Therefore Fi^^ = Fq^i and by connectedness of t we then get Fi^ = Fqx for all x in A, 
i.e., di = (Iq = d2- It then easily follows from (3e.7) that also ni = uq = n2- 

We are left with the case where the off-diagonal part t' is bipartite. In this case we get from 
(3e.9): 

\6oix)\-\diix) - doix)) = ci-ir (3e.l0) 

for some constant c. We shall show now that if c 7^ then Dq cannot be minimizing. By 
connectedness of t we see that (3e.8) implies that 

6i{x) = ai6o{x), (3e.ll) 

where oi is a constant (recall that here ^0(2;) is positive) . 
We turn to the diagonal part of (3e.7): 

(t^x -11- Uxd^{x)f + Ul\5r{x)\^ = {t^x -II- U^dof + Ul5^{xf. (3e.l2) 

If we insert (3e.l0) and (3e.ll) into (3e.l2) we get 

2c{-lYU-^5o{x)-\txx -II- U^doix)) + + |ai|2 = 1. 

We therefore conclude that if c 7^ then 

ixx - M - U^doix) = {-irUJoix)C, (3e.l3) 

where C is a constant. We can restate this as 

/t' \ /(-irc 1 \ 

T,-UDo=(^ _ J - C/.<5o(.) ^ -(-ircj- ^'"''^ 
We can choose a pseudo-spin rotation W as in (3d.52) such that 

Wi = { ^ . (3e.l5) 

Since t' is assumed to be bipartite we see that the unitary operator 
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has the effect VW{T^ - UD)WW^ = W{T^ - UD)Wl 

From (3b. 31) and (3c. 4) we see that if V is the minimizer corresponding to Dq then T' = 
F^{{T^ - UDq)) where Fp is the function F^(A) = [1 + exp(/3A/2)]-i (for /? = 00 we have F^{X) = 
hm^^oo Fp{X) = x(A)). Hence VWT'WW^ = WV'W^. For the potential Do this imphes that 



Since WDqW"^ is traceless we see in particular that the diagonal entries of WDqW^ must be zero. 
If we write 



we see by comparison with (3e.l5) that wc must have dQ{x)5o{x)~^ = ( — 1)^C It then follows 
from (3e.l3) that t^x — ^ = 0. Since we are in the case where t' is bipartite this implies that t — /i 
is bipartite contrary to our assumption in the lemma. ■ 
We can now state the main result of Chapter 3. 

3.12 MAIN THEOREM (Characterization of HF states): The HF Gibbs states or 
ground states are unique modulo global gauge transformations. More precisely, all HF Gibbs states 
at the same inverse temperature /3 or all HF ground states (i.e., the minimizers of the functional 
(3b. 23) and (3c. 2) respectively) are related by global gauge transformations. 

Proof: By Theorems 3.3 and 3.4 there is a one to one correspondence between the HF states 
described by T' minimizing —Vp,^ and the potentials D minimizing 7^/3,^- The uniqueness of D 
modulo gauge transformations follows immediately from the results about "7^/3,^; Lemmas 3.6-3.8 
and Lemmas 3.10 and 3.11. I 

We remark that Theorem 3.12 allows us to characterize the HF Gibbs states by the global 
gauge and the inverse temperature. Indeed, the correspondence between (3 and the HF Gibbs state 
T'{(5,w) is continuous in [3 for every fixed global gauge w (w being a pseudo-spin rotation, or a 
phase change). In particular, this mapping is continuous at /3 = 00 and, a posteriori, we could have 
spared the entire discussion of the zero temperature state by simply appealing to this continuity. 

We shall see in Sect. Ill.g that if the temperature is high enough their maybe no symmetry 
breaking and hence the Gibbs state is unique. There is no degeneracy caused by gauge transfor- 
mations. If however, the temperature is small it follows from Lemma 3.5 that there is symmetry 
breaking for any finite system. In the translation invariant case studied in Sect. Hl.g, the symmetry 
breaking phase transition persists in the thermodynamic limit. 

Ill.f SPATIAL SYMMETRIES 

In Sect. Ill.d we studied global gauge symmetries of the Hamiltonian. In the present section we 
shall address the question of spatial symmetries, i.e., symmetries of the lattice or, more precisely, 
of the t matrix and the coupling constants Ux- 

As a special example, in the next section we shall explicitly determine the minimizing matrix- 
valued potential D when the lattice is translation invariant and the coupling constants are 
independent of x. 
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By a spatial symmetry we understand an invertible transformation r : A hh' A of the lattice 
A. We say that t is magnetically invariant under r if tliere exists a map a-r A. <—>■ H such that 

tT{x)T{y) = ex.p[-i{ar{x) - ar{y))]txy 

Equivalently, if we reahze t as an operator on H magnetic invariance means that t commutes with 
the unitary transformation m(r) on H defined by 

m{T)\x,a) = exp{—iaT-{x))\T{x),a). 

We need the phase factor a in order to treat non-real t. If t is real we may of course choose 
ar = 0. As an example consider a two-dimensional torus, i.e., a finite cube A in with periodic 
boundary conditions. The original (real) hopping matrix considered by Hubbard, Kanamori and 
Gutzwiller is invariant under pure translations (m(T) with = 0) on the torus. Consider, however, 
the (complex) hopping matrix t which differs from the elements of the original matrix by the 
multiplication of complex phases and which correspond to having a fixed magnetic flux through 
each unit square. Then t is not invariant under pure translations but rather under the magnetic 
translations m(r) which are compositions of translations and gauge transformations {a^ 7^ 0). 

The family of all transformations r for which t is magnetically invariant and the coupling 
constants satisfy {/^(a,) = Ux for all x in A naturally forms a group Q which we call the spatial 
symmetry group (of t and U). Notice that r 1-^ m{T) need not be a unitary representation of 
Q. In fact for the two-dimensional torus the translations commute while the magnetic translations 
do not. (If T\,T2 denote the translations of unit length along the first and second direction respec- 
tively we have m(Ti)m(r2) = exp(iij!')m(r2)?7?.(ri), where (j) is the flux through the unit squares.) 
We emphasize that the group Q itself need not be abelian, it could be one of the crystallographic 
groups (e.g. if we study the Hubbard model on the lattice formed by the carbon atoms in the 
Buckminsterfullerene Ceo molecule). We say that t is translation invariant if the spatial sym- 
metry group Q acts transitively on the lattice A, i.e., if for any two points x,y & A there is a r in 
Q such that t{x) = y. If t is bipartite (and connected) it is easy to see that each clement r of ^ 
must either map the A and B sublattices into themselves {t(A) = A and t(B) = B) or map the A 
sublattice to the B sublattice {t{A) = B). It is clear that the latter type of transformations exist 
only if \A\ = \B\. 

Corresponding to r G ^ we define a Bogoliubov transformation Wr on the Fock space T by 

Wrcl^Wl = exp(-ia,(x))4(^)^^. 
The unitary matrix Wr corresponding to Wr is 

V rn{T) J 

If r belongs to the spatial symmetry group then the Hubbard Hamiltonian is invariant under 
Wr- We shall prove in the next lemma that if t is also non-bipartite then the HF Gibbs states 
are also invariant under Wr- In the bipartite case the states are also invariant under transfor- 
mations such that t{A) = A and hence t{B) = B, they need, however, not be invariant under 
transformations such that t{A) = B. 
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3.13 THEOREM (Spatial invariance): Let Q he the group of spatial symmetries oft and 

U. 

(a) If t is bipartite (but not necessarily real) and fi = all minimizing potentials D for TZp^ 
satisfy that n(x) = ^Tr[D'^] is invariant under Q, i.e., n{T{x)) = n{x) for all t G Q. All 
minimizing T' for —Vf^fi are invariant under unless r maps the A sublattice to the B sublattice 
(t{A) = B). As a consequence the HF ground and Gibbs states are invariant under Wr unless 
t{A) = B. 

(h) Assum,e t real and t — fi not bipartite Then the minimizers D for T^js,^, T' for —V/s^n and 
the HF ground and Gibbs states are invariant under Q . 

Proof: Part (a) follows from the strict convexity proved in Lemma 3.10 since it implies that n 
is unique and hence invariant. The second statement in part (a) is a consequence of Lemma 3.7. 
where all the possible minimizers D are described. In fact, the minimizers D are invariant under 
r unless t[A) = B because the 2 x 2-matrix valued function Wx defined in Lemma 3.7 is constant 
on the A sublattice and on the B sublattice. 

Part (b) follows from Lemma 3.11 and Lemma 3.6 ■ 

Ill.g THE TRANSLATION INVARIANT CASE 

In this section we shall explicitly determine the minimizing matrix-valued potential D under 
the additional assumption that the lattice is translational invariant, i.e., that the spatial isometry 
group G of t and U acts transitively. In particular this means that U is constant independent of x. 

As shall be shown, there occurs a phase transition as the temperature varies about a critical 
value Tc, or /?c, respectively, provided 1 < |^Tr-^^[|t — /^|~^] < 00 holds and assuming that the 
spectrum of t — is symmetric about 0. More precisely, we will show that = 0, i.e. n{x) = 
for all X G A, if /3 < /?c and n{x) equal to a non-zero constant no for all x G A in case of /3 > f3c- 

Since, as we pointed out in section Ill.d, n{x) = no > goes along with having non-vanishing 
aligned pseudo-spins everywhere in A, this implies long range order which is off-diagonal in case 
that 8{x) ^ for ah x G A. 

Secondly, due to the analytical dependence of the pressure on Z), the transition from n = for 
/3 < /3c to n > for /3 > /3c also indicates that the pressure is nonanalytic at /3c- Notice that in our 
model there is no thermodynamic limit |A| — > 00 required to yield a nonanalytic thermodynamic 
potential. 

3.14 THEOREM (Pressure in translation invariant case): If t (with as usual t — 
either bipartite or real) and U are translation invariant the pressure is 

V{I3, n) = - min I -2/3-^Tr [in cosh | ^{t - nf - 2Ud(t - fi) + U^v] +Uri\A\} 

Cp<.ri L L J J 

+ (2/3-Mn2 + /x)|A|, (35.1) 

where the minimum is over real constants d and rj. The minimium in (3g.l) occurs at unique 
values do {P, fj.) andr]o{P,n) satisfying 

do{P,fif <Vo{P,fi)<i (35.2) 
Ift — jj, is bipartite then do(/3,/x) = 0. 
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Moreover, if D is any minimizer for 1^/3, n then for all x we have n{x)^ = ^Tr[D'^] = rj. If 
t — ji is not bipartite we also have that the upper diagonal element of is d{x) = do for all x. 

Proof: Since a; t— > — In cosh ^/x is strictly convex and since we are minimizing over the convex 
domain {dg < ryo} we conclude that the minimum occurs at unique values. This is just a special 
case of the more general statements in Lemmas 3.10 and 3.11. 

If t — /X is bipartite it is unitarily equivalent to — (t — //) it then follows from (3.113) that if 
do and ?7o are minimizing values then so are —do and 770. By uniqueness we therefore have that 
do = 0. 

If t — /X is bipartite we may without loss of generality assume t bipartite and /x = 0. It then 
follows from part (a) of Theorem 3.13 that all minimizers D of 7^;3,^ satisfy that n{xf = \Tt[DI] 
is a constant rj. We see from Lemma 3.7 that is of a form such that (T^ — UD)^ = + U^D"^. 
Since the minimum in (3g.l) occurs for do = we see from (3b. 30) that (3g.l) is indeed a correct 
formula for the pressure. 

If i — is not bipartite but real and if the functions d and rj define a minimizer for TZis^n (in 
the sense of Lemma 3.11) it follows from the uniqueness proved in Theorem 3.5bb that d and 77 
are independent of x. It is then clear from (3b. 30) that (3g.l) is correct and that the parameters 
r] = r]o and d = do are the unique minimizers in (3g.l). ■ 

We now restrict to the case when f — is bipartite or more generally to the case where t — fi 
is unitarily equivalent to —{t — n) and 

1< ^TV[|t-/x|-i] <oo (35.3) 

holds. Then a critical inverse temperature < /3c < 00 is uniquely determined by 



1 = — Tr 
|A| 



li\ ^ tanh^/3c|t — 



(35.4) 



since the right side in (3g.4) is continuous and grows monotonically with f3c from to U\k\ ^Tir^t — 

The following theorem establishes that 13c is, indeed, critical, provided the spectrum of t — /j, 

is symmetric about because, then, do{(3,^) always vanishes and the issue of determining the 

minimizer D simplifies. Notice that in the case where do = 0, the gap around zero in the spectrum 

1/2 1/2 
oiT^ — UD is at least 2Ur]Q . If t — /x has a zero eigenvalue then 2UrjQ is precisely the value of 

gap. 

3.15 THEOREM (Gap equation): Let t be real, translation invariant and Ux = U > for 

all X E A. Assume that t — n and ji — t are unitarily equivalent as in the case of a bipartite lattice 
with 11 = 0. Let (3c he given by (3g.4) in case (3g.3) holds and (3c 00 otherwise. Then % = 
for (3 < (3c. Moreover, r]o{(3,fi) is a strict monotonically increasing function in (3c < (3 given by the 
gap equation 



ii^Tr 
' |A| 



[{t- fif + U^r,o) '^%anh(^fV(^-/^)' + f^'%) 



(35.5) 



provided (3c < 00. 
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Proof: Equation (3g.5) follows by setting the derivative with respect to 77 of the expression on 
the right side of (3g.l) equal to zero. Note that since tanhx < 1 (3g.5) implies 

in agreement with our previous condition ijq < ^. 

The monotonicity of r)o{P, //) in /3 > (3^ is straightforward from (3g.5). ■ 
On a translation invariant lattice we may represent t — /x by its eigenvalues £k for k G BZ, 
the Brillouin zone, whose volume we denote by \BZ\. In the thermodynamic limit, |A| — > 00, the 
gap equation (3g.5) then takes the more familiar form 

1 = J^^{el + U^m)''^' tanh(^f ^47c^)dk, 

which is the BCS gap equation. 



1< i^TV 
- 2 lAI 



IV. THE GENERALIZED HF THEORY FOR THE 
HUBBARD MODEL WITH REPULSIVE INTERACTION 

IV.a LINEARIZATION OF THE PRESSURE FUNCTIONAL 

In this section the generalized HF theory will be applied to the Hubbard model with repulsive 
interaction. We continue to use the notation of Chapter III and consider the Hamiltonian 

x,y£X xGi 
a 

which differs from iJ_ in (3a. 1) in the reversed sign of the interaction, i.e., we again assume Ux > 0. 
A close look at the pair interaction reveals that 4 I'^x t4 i^x l ~ t4 I'^x i^x f ^ 0- Thus the 
interaction (corresponding to the operator V in Chapter II) is repulsive and Theorem 2.11 applies. 

Hence, we may restrict our attention to 1-pdm of the form T = ^ _^ and the energy 

expectation reduces to 

^(r) = TrK7] + E Ux{b^{x) - i][7i(x) - i] - \Mx)f}, (4a.2) 
where 70- (a;) := {x,a\'y\x,a), 7*(a;) := {x,1 |7|x,|) and we denoted 

n-Co" (^-^) 

on H. Obviously, £{T) depends on 7 only and we will write £"(7) := S{r). In fact, dealing merely 
with 7 we will have to consider only operators on H rather than H^BH throughout this chapter. 

We denote ■= IxQ^x for any operator Q on H where Ix is now the projection onto 
functions in H vanishing everywhere except at x. 
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Equipped with this notation we may write 



7x = Ix (4a.4) 



and one easily verifies that 

[7t(^) - i][7i(^) -W- Mx)? = \ (Tr[7. - ^1.])' - ^Tr [(7. - \r.f\ . (4a.5) 
The entropy depends merely on 7, too, namely 

5(7) = -5Tr[r inr + (1 - r) ln(l - r)] 

= -Tr[7ln7 + (l-7)ln(l-7)], (4a.6) 
and so does the pressure expectation 

-Vfi,^{^) =IV[r;7] + I E (^[^-] - 1)' - i E ^-Tr [(7x - 

+ /?-^TV[7ln7 + (l-7)ln(l-7)], (4a.7) 

where we again denoted £'(7) := <S'(r) and V^^^i^i) := Vjs^^iiT). For the examination of the HF 
ground states we denote Voo,n lim/3^oo 'P/3,ij. = £ — jJ-N. Let us introduce an auxiliary functional 

-VpAl) =Tr[r;7] - \ E ^-Tr [(7x - \p{x)l.f] 

X 

+ /3-^TV[7ln7 + (l-7)ln(l-7)], (4a.8) 

denoting := Tr[72:]. Of course, we write Voo,n ■~ lini/3^oo "P/J.p- Notice the formal similarity 
between (4a.8) and (3b. 23). Replacing by and V by 7 in (3b. 23) we arrive at (4a.8) excpet 
that p{x) is missing (and that there is an extra unimportant term — At|A| in (3b. 23)). The reason 
p{x) is missing in (3b. 23) is that Tr[r^] = 1 is a consequence of the form (3b. 13). This is the 
crucial formal difference bewteen the attractive and repulsive cases. We claim that 

-P/3,m(7) = -r^Ai) + -4 E ^M^) - 

X 

Indeed, since 'jx — ^p{x)lx has zero trace, 

-VpAl) + "PpAl) = k E U.{Hx) -1]'-Tr [(7x - ^l.f] + Tr [{j, - ^p{x)lxf] } 

X 

= -2 E M^P(^^ - 1]' - TV [i[p(x) - i]Hx] } 

X 

l^UM^)-!]''. (4a.l0) 



4 . 

X 



Equation (4a.9) is an important observation. With our machinery developed in Chapter III 
we can only determine the minimizers for —Vp^^ rather than —Vp,^. This substitution is justified 
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only if we can show that a 1-pdm 7 minimizing — P/3.^t also minimizes —Vp^^. Wc succeed in doing 
so only in case of bipartite hopping matrices t with chemical potential /Lt = 0. In fact, in this case 
we will prove that 

p{x) = Tr[7,] = C7|7|x, a) = 1 (4a.ll) 

for all lattice points x G A. We shall refer to (4a. 11) as the constant density lemma because 
of its similarity to the main theorem in [LLM]. Equation (4a.ll) establishes the formal analogy 
between the attractive and repulsive cases. 

We start with the analysis of — and its minimizer. We will denote 

-V{f5,^l)■.= min -Vp,^{^) (4a.l2) 

0<7<1 

for positive or infinite /3. In analogy with (3b. 25) we first observe that 

-Tr [(7^ - \p{x)l^f] = min {-2IV[L>^7,] + Tv[DI]] , (4a.l3) 

a,o 

/d{x) 6{x) \ 

where := . Because of (4a. 13) the proof of the following lemma is a line-by-line 

\6{x) -d{x)J 

copy of the one for Theorems 3.3 and 3.4. 

4.1 LEMMA: For all < P < 00 and all fj, we can write the auxiliary functional P{(3,fi) as 
the following variation over the functions d and S. 

-r{p,n) = mmTZf,,^{d,6) -2p-'\A\\n2, (4a.l4) 

d,o 

where 

np,^{d,6) ■.= np,^{D) :=-/3-iTV[ln cosh f(T;-C/Z))] + |J^[/,Tr [1)2], (4a.l5a) 

X 

and 

n^Ad, S) := n^AD) ■= - iTr|r; -UD\ + ^Y1 ^-^1^1]. {4.al5b) 

X 

If a potential D minimizes 'R-js^^ then the operator 

(l + exp /?(r; -UD)^ , forp<^ or 7 = xi^ - UD), for p = ^ (4a.l6) 

minimizes the auxiliary functional —Pp,^ and satisfies the consistency equation 

7x = {D^ + ^Tr[7,]) 1,. (4a.l7) 

Conversely, if j is a minimizer for —Pp^n then the potential D defined by (4a. 17) minimizes Ti-p^^i 

and satisfies (4a. 16). 

As in the attractive case the minimum in (4a. 14) occurs for functions d and 5 satisfying 
d{x)^ + \5{x)\'^ < I for allxe A. 

We prove an analog of Lemma 3.5 by merely replacing by T'^ in the proof of that lemma. 
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4.2 LEMMA (Gap estimate): Let D be a minimizing matrix-valued potential for TZf^^fj, or 
T^oo,ii and denote the eigenvalues ofT^ — UD by ei, 62, • • • , e2|A| • Then, for any j = 1,2, . . . , 2\A\, 

we have 

|e,| > iC/„,i„|A|-i -/3-i21n2 or |e,| > it/^i„|A|-\ (4a.l8) 

respectively, where Umin '■= niinlt/^;}. 

Lemma 4.1 allows us to concentrate now on the determination of the functions d and 6 that 
yield a minimizer D for 7^/3,^ or T^oo,^ respectively. Note that TZis^n in the repulsive case diflFers 
from TZ/3^n in the attractive case only by the replacement of by T^. 

IV.b CONSTANT DENSITY LEMMA FOR BIPARTITE LATTICES 
AT HALF-FILLING 

In this section we will assume that the hopping matrix t is bipartite (possibly non-real) and 
the chemical potential fi equals zero. It will turn out that this choice of t and /i allows us to 
conclude that Tr[7] = jA| = ^dimH which is the reason for calling this case half- filling. 

4.3 LEMMA (Antiferromagnetic spin alignment): Let t be bipartite (but not necessarily 
real) and /x = 0. Then 

n^Ad,^) > n^,^ ((-i)Vd2 + |,5|2,o) . (46.i) 



Lf n{x) := ■yjd'^ix) + |(^(a;)P > for all x & A equality holds in (4b. 1) if and only if for all x ^ X 

( d{x) d{x) \ 



S{x) —d{x) 



{-lYn{x)w w'^ (46.2) 



for some unitary 2 x 2-matrix w independent o/x G A. 

The proof of Lemma 4.3 is in complete anology to the proof of Lemma 3.7. We merely have 
to replace Tq by Tq. Notice that (4b. 2) results from 



djx) _ d{y) 
n{x) n{y) 



O5 txy 



Sjx) _ S{y) 
n{x) n{y) 



0, (46.3) 



which replaces the conditions (3d. 25a) and (3d. 25b) in the proof of Lemma 3.7. Note that (4b. 2) 
shows that the potential D has the staggered order characteristic of antiferromagnetism. The form 
(4b. 2) is exactly what we need to prove the constant density lemma. 

4.4 LEMMA (Constant density): Let t be bipartite (but not necessarily real) and n = 0. 
Define 7 by 

7 := (1 + exp[/?(T^ - UD)]y\ (46.4) 
where is subject to condition (4b. 2) for some unitary 2 x 2-matrix w. Then 

Tr[^x]=Y.{x,a\^\x,a) = 1. (46.5) 

In particular, Tr['y] = J2x "^[ix] = |A|. 
Proof: We want to show that 

= J^(x,c7|i -7|x,c7) = U](x,c7|tanhi/3(r^ - C/L»)|x,(7). (46.6) 
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We define a unitary transformation 

V^:=w(^l -^^wH-iy. (46.7) 

One easily checks tliat 

- UD)Vl = - UD). (46.8) 

Thus, 

^{x,a\i&ah\p{TQ-UD)\x,a) = ^{x,a\V^ [tanh i/3(ro - C/L>)] Vl\x,a) 

u a 

= - ^{x, (t\ tanh i/?(T^ - UD)\x, a), (46.9) 

cr 

which is equivalent to (4b. 6). ■ 

Now we are in a position to determine the actual pressure. 

4.5 MAIN THEOREM (Ground state and positive temperature pressure): Let t 

he bipartite and fi = 0. Then for all < P < 00 

--pi^p^O) = -r{p,0) = mmTZ0fi{d,S) -2l3-'^\A\\n2. (46.10) 

d,5 ' 

If a potential D minimizes 7^/3, for < f3 < 00 then 7 defined in (4a. 16) is a HF Gibbs state 
(P < 00) or ground state (P = 00). This 7 satisfies the consistency equation 

7x = {Da, + i)lx. (46.11) 

The function n{x) = Y^ci(xp + |^(x)P corresponding to a minimizer D is unique and vanishes 
either everywhere or else nowhere. The potential Dx must he of the form (4b -2) for some unitary 
matrix w. As a consequence, the HF Gibbs states and ground states are unique modulo gauge 
trans fo rm ations. 

Proof: The right side of (4b. 10) is clearly a lower bound to —V{P, 0) since — P/3,o(7) ^ — ^/?,o(7) 
for any < 7 < 1. On the other hand by Lemma 4.3 the minimum on the right side of (4b. 10) can 
be attained by a potential D satisfying (4b. 2) (e.g. with w equal to the identity matrix). We then 
have a minimizer 7 for V defined in terms of D by (4a. 16) (with = 0). Note that because of the 
gap estimate we know that x(^o ~ UD) = lim^^oo(l + exp[/3(TQ — UD)])~^. It therefore follows 
for both finite and infinite P that 7 satisfies the constant density relation (4b. 5). Thus by (4a.9), 

-V{p, 0) < -Vp,o{l) = -P/^Al)- (4b.l2) 

This conludes the proof that 7 is a HF Gibbs or ground state. The remaining part of the heorem 
follows by an analysis identical to the one leading to Lemma 3.8, Lemma 3.10 and Theorem 3.12 
(part of the conclusion is that 7?./3,o((— l)^?^; 0) is convex as a function ofrj = and strictly convex 
at its minimum). ■ 

Finally, we mention that a result on spatial symmetry analogous to Theorem 3.13 (with the 
same proof) holds in the repulsive case. 

4.6 THEOREM (Spatial invariance): Assume t is bipartite (not necessarily real) and 
fi = and let Q be the group of spatial symmetries oft and U. For any minimizer {d, S) ofTZp^o we 
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have that the function n{x) = ^/ d{xy + \S{x)\'^ is invariant under Q. The HF ground and Gibbs 
states are invariant under Wr for t & Q unless t{A) = B. 

IV.c PARTICLE-HOLE SYMMETRY 

In the special case of a real bipartite hopping matrix t at half-filhng /x = there is a more 
elegant way of deducing Theorems 4.5 and 4.6 from Theorems 3.2, 3.5 and 3.8 by means of a partial 
particle-hole transformation only on the j-spins. More precisely, we define the Bogoliubov 
transformation Wp/i,t by 

Then 

x,y X 

~ y^A^xy^xA^vA + ^^y^x,l^y,0 ~ 5^ 

x,y X 

X 

=H_, (4c.2) 

using txy = txy = tyx, the canonical anticommutation relations and ^^txx = 0. Since the 
Bogoliubov transformation leave the set of quasi-free states invariant, it immediately follows that 

E^^{H_) =mi{p{H_) I p is quasi-free } 

= inf {pwp^^(if_) I p is quasi-free } 
= inf {p{H+) I p is quasi-free } 

=E^^{H+) (4c.3) 



and a similar equality holds for the pressure at finite temperatures. Thus, there is clearly a one- 
to-one correspondence between the l-pdm r_|_ minimizing the generalized HF pressure functional 
for and the 1-pdm r_ minimizing the generalized HF pressure functional for H- via 

Note that if t is not real the attractive and repulsive cases are not unitarily equivalent. In 
fact, it follows from Theorems 3.5 and 4.5 that in the attractive case there may be no more than 
two minimizers, while in the repulsive case there is a continuous family of minimizers related by 
spin 517(2) transformations. 



IV.d FERROMAGNETISM AT INFINITE REPULSION 

In the preceding sections we found that for bipartite lattices at half-filling, N = |A|, the 
minimizing 1-pdm always has antiferromagnetic order, i.e., the sign of d on the A-sublattice is 
always opposite to the sign on the B-sublattice. On the other hand Nagaoka's Theorem states 
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that the true ground state (which happens to be the free particle HF state) for A'^ = |A| — 1 has 
maximal spin S = N/2, provided one takes U = 00 (see [NY], [TD], [TH]). In our language this 
would mean that the ground state 7 would satisfy (after an SU{2) rotation) 7j,(x) = 0. Setting 
U = 00 mathematically, means, projecting out the vectors with doubly occupied sites from the full 
Fock space. 

In the context of the generalized HF approximation one can also make sense of C/ = 00 and 
we will derive an analog of Nagaoka's Theorem for the HF minimizer. The analogy goes too far, 
however, because as we will prove for any 1 < A'' < |A| — 1 the HF ground state 70 will be fully 
spin polarized. In contrast, for the true Hubbard model this does not hold, i.e. for N = \A\ — 2 
(see [DW], [DFR], [TB],[SA]) or for < 0.51|A| [SKA], the ground state does not have S = N/2. 

Let us start by defining what is meant by infinite repulsion in the context of generalized HF 
theory. Recall from (4.2) that the energy functional for positive coupling Ux = U > becomes 



Sij) =TV[ra] + [/ J] {[7^(x) - i][7i(x) - 1] - |7*(x)|2} 
xeA 

We consider the particle number N = Tt[^] fixed. The last two terms in (4d.l) are therefore 
constants that we may ignore when determining the minimizing 7. The limit U ^ 00 yields the 
constraint 

7T(.T)7i(x)-|7*(x)|2 = (4d.2) 
for any x G A. More precisely, if we define the Hartree-Fock energy of N electrons by 

E^^iN) := inf{^(7) + ^UN - iC/|A| I < 7 < 1, Tr[7] = N^ (4d.3) 

then limc/^oo E^^{N) = E^l^{N) where 

E^l^{N) := inf|Tr[T^7] < 7 < 1, Tr[7] = N,j fulfills (4d.2) }. (4d.4) 

We remark that N < |A| is automatic in (4d.4) because the constraint (4d.2) is equivalent to 
Det[7a;] = 0, which together with < 7 < 1 implies that Tr[73,] < 1 and hence N = Tr[73;] < 
|A|. 

4.7 THEOREM (Ferromagnetism at infinite U): Let ei < 62 < . • • < e\\\ denote the 
eigenvalues oft. Then 

N 

E^l^{N) = Y,e^. (4d.5) 

1=1 

A family of minimizing 7 for the variation in (4d.4) corresponds to the ferromagnetically saturated 
states i.e., is given by 

7 = w;(^^ ^^wl (4d.6) 

where w is any SU{2) matrix and Pn is the spectral projection onto the eigenvectors with eigen- 
values ei, e2, . . . , Civ. 
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N 



Proof: It is clear that E^^^{N) < ^ e^. Because any matrix of the form (4d.6) is admissible 

i=l 

for the variation in (4d.4). On the other hand given any 7 satisfying (4d.2), < 7 < 1 and 
Tr[7] = N. We can write 



7 



^\'Pi){'Pi\ = ^\fi®9i){fi®9i\, 



(4d.7) 



where {ipi\ipj) < 5ij for all On the right side of (4d.7) we regarded H as (B Hi, and we 
denoted fi{x) := (pi{x, j) and gi{x) := ipi{x, [). Now, define on Ha 



■=Y^\fi){fi\ + \9i){9il 



(4d.8) 



such that 7 := 



7t 







is completely spin-up polarized. The l-pdm 7 naturally fulfills (4d.2). 



Our goal is to show that 7| < 1 and hence 7 < 1 since this, together with < 7 and Tr[7] = A'^ 
implies (4d.5) because 



JV 



TVK7] = TV K7] = Tr[tj^] > J2 ^i- 



(4d.9) 



In order to demonstrate 7| < 1 we need to take a closer look at the constraint (4d.2) first. 
We may rewrite (4d.2) as 



(4d.l0) 



for all X G A. Let us denote 



Ao := G A I fi{x) = f2{x) 



/iv(x)=0}. 



(4d.ll) 



On Aq (4d.l0) holds trivially, but on the complement it yields the existence of a complex number 
a{x) such that 

gi{x) = a{x)fi{x), g2{x) = a{x)f2{x), ...,gN{x) = a{x)fN{x) (4o?.12) 
for any x G A \ Aq by Schwarz' "equality" . We define a normal operator A, AA^ = A^A on Ha by 

A:= a{x)\x){xl (4d.l3) 

xeA\Ao 

and denote the projection onto Aq by B := YlxeAo Hence, we may express (4d.l2) as 

gi = Afi + Bgi. By means of A and B we can write 



S^J > {^^\^J) = iMf,) + = + A^A)fj) + {g,\Bg,) 

={m + AA^)f,) + {g,\Bg,). 



(4d.l4) 
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We define Fi := \fi){fi\, Gi := \gi){gi\, Ri := B\gi){fi\. Now, 7t < 1 and, therefore, (4d.5) are 
direct consequences of the inequality 

7T -7? > {y.[A,F,\ + r}i {y.[A,F,\^ + r\^ + (Y.^\^ (E^^) • (4^-15) 

The inequality (4d.l5) follows from (4d.l4), 7^ = Z^j -^i + Gj, and AA'^ = A since 

FiF, =\h){h\f,){f,\ < I/.) (5,, - iMA^Alf^) - {g,Bg^)){f^\ 
=SijFi - FiA^AFj - R\Rj,= 5,jFi - FiAA^^ F^ - rIRj, 

GiGj =\9i){9i\9j){9j\ < \9i){Sij - {fi\fj)^{gj\ 

=SijGi - AFiFjA^ - RiFjA^ - AFiR] - RiR], (4d.l6) 

FiGj =FiAFjA^ + FiARt . ■ 

The previous theorem does not tell us when the ferromagnetically saturated states are the only 
HF ground states. The proof only gives the following criterion, which, unfortunately may be very 
difficult to verify. Assume that < Gn+i which ensures the uniqueness of Pn- The criterion for 
uniqueness of the ferromagnetically states given by 1-pdm of the form (4d.6) is that the spectral 
projection Pjv be connected in the sense that any two points x,y £ A can be connected by a path 
X = xq, xi, . . . ,Xn = y for which (xj|P/v|2;j+i) 7^ 0. To prove this criterion suppose that 7 in (4d.4) 
is a minimizer, i.e., 

TV 
i=l 

where 71 is defined as in (4d.7). Hence, 71 = Pn = 7| and (4d.l5) implies 



J^i?^ = 0, 



(4d.l8) 



which, in turn, gives us 

^JV = 7T = E + ^» = E + ^^^^^ + + ^1^^ + ^^i^ 

i i 

= ^Fi + AFiA'' + BGiB. (4d.l9) 

i 

The right side of (4d.l9) is not connected between Aq and A \ Aq unless Aq = or Aq = A. The 
case Aq = A simply means that 7 is completely spin-down polarized, in accordance with (4d.6). 
Conversely assuming Aq = 0, we observe that (4d.l9) also implies [A, Pjv] = and obtain 

[a{x) - a{y)] {x\PN\y) = (4d.20) 

for all X, ?/ G A \ Aq = A. But Pjv was assumed to be connected and, therefore, a{x) =const for all 
X e A. In other words, fi = agi for all i which implies (4d.6). 
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V. SUMMARY OF HF THEORY OF THE HUBBARD MODEL 

V.a INTRODUCTION 

Our aim here is to give some perspective to the results in Sects. Ill and IV by summarizing 
them — with special emphasis on symmetries and their breaking. The first task is to define sym- 
metry breaking. We start with a Hamiltonian H that, in many cases, is invariant under some 
symmetry group, G, each element of which is represented by a unitary operator on our Hilbert 
(Fock) space of dimension 41^^!. The representation of the group G might be a ray representation, 
as in the case of "magnetic translations" . The unitary operator corresponding to an element w e G 
will be denoted by W{w) on the one-particle space H (BH and by W{w) on the Fock space J^. 
Since the only unitaries that transform HF states into HF states are Bogoliubov transformations, 
we restrict our attention to groups consisting of such transformations. This restriction is not really 
a cause for disappointment because all the symmetry groups that are usually considered, such 
as rotations in real space, rotations in spin space, translations, etc. are, in fact, represented by 
Bogoliubov unitaries. The reason is simply that most symmetry groups in physics are defined by 
their action on one particle and are then extended to iV > 1 particles by tensoring. This is exactly 
what number conserving Bogoliubov transformations do. 

We shall take the point of view that a state breaks a symmetry if the state is not invariant 
under the action of the corresponding symmetry group. This is formalized as follows. 

5.1 DEFINITION (Broken or unbroken symmetry): Let G be a group represented by 
unitaries as described above. Let p be a state on the operators on Fock space T corresponding to 
the Hilbert space 7i. We say that the G symmetry is unbroken in the state p if, for each w G G, 
p is invariant under the Bogoliubov transformation W{w), i.e., p{W{vj)AW{wy) = p{A). If the 
state is not invariant for allu) e G we say that the G symmetry is broken in the state p. 

We are particularly interested in HF ground states or finite temperature Gibbs states for a 
Hamiltonian H invariant under the action of some group G, i.e., W{w)HW{wY = H for all w & G. 
The interesting question is then whether or not the G symmetry is broken in these states. 

The contrast between the usual theory and the HF theory should be kept in mind. While the 
original Schrodinger equation Hi/j = Eip defines a linear theory, the HF theory that approximates 
it is intrinsically a nonlinear, one-particle theory. Linear combinations of HF wave functions 
(which are Bogoliubov transforms of the zero-particle vacuum) are not HF wave functions. On 
the level of states (either pure states, A 1-^ {iIj\A\iIj), or Gibbs states), the HF states do not form 
convex sets. Usually, if pi and p2 are states (either ground states of some Hamiltonian H or Gibbs 
states of H at some temperature T — in the thermodynamic limit there can be more than one) then 
p = Api + {1 — X)p2 is an admissible state (ground or Gibbs). As discussed in the introduction 
this is not true for HF states because p is not usually a HF state when pi and p2 are. In the usual 
theory we can ask for the extremal states (i.e. those states that are not convex combinations of 
other states at the same temperature) and ask about their properties with respect to symmetry 
operations. An example to keep in mind here is the Heisenberg Hamiltonian, for which there is a 
magnetized ground state. By taking a convex combination of all ground states one can construct 
a ground state that is SU{2) invariant, but this state is not extremal in the set of ground states. 
Prom this example we learn that in the usual theory symmetry breaking should be sought only 
with extremal states — which correspond to pure phases. 

The HF states, on the other hand, do not form convex sets and therefore we cannot talk 
about extremal states. We regard each HF state, heuristically, as playing the role of an extremal 
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state. Indeed, in the usual theory symmetry breaking in a finite system is infrequent; typicahy it is 
necessary to pass to the thermodynamic limit in order to see it. In contrast, if symmetry breaking 
occurs in HF theory, it is usually manifest for the finite system. The following discussion refers to 
either the ground state or to positive temperature states. It is to be understood that the symmetry 
breaking displayed in our three tables may not actually occur. In particular, they will usually not 
occur if the temperature is high enough. 

V.b SYMMETRIES OF THE HUBBARD HAMILTONIAN 

To begin our summary of symmetry breaking in the HF theory of the Hubbard model we first 
list the symmetries of the Hubbard Hamiltonian. There are two types of possible symmetries of 
global gauge symmetries and spatial symmetries. The spatial symmetries of H depend on the 
spatial symmetries of t and U as explained in Sect. Ill.f. The gauge symmetries depend on the 
presence or absence of the following three properties: Bipartiteness of t, reality of and /U = 0. 
Whether the gauge symmetries are broken may depend on the sign of the interaction. 

(i) Spatial symmetries: In Section Ill.f, Theorem 3.13, we found that the spatial symmetries 
(if any are present) are broken only for a bipartite t and then only when /i = 0. Each spatial 
symmetry either maps A to A and i? to -B or maps Aio B and B to A. Even for bipartite t and 

= 0, the A to ^4 symmetries are never broken. It is only the transformations r in the spatial 
symmetry group which take the A sublattice into the B sublattice, i.e., t{A) = B (hence of course 
t{B) = A) that are broken. In fact, in this case we must require that the two sublattices A and 
B have the same number of points {\A\ = \B\). Whenever we refer to a broken spatial symmetry 
in the tables below, we always mean that it is only the symmetry of maps from A to B that is 
broken. (Note that there can be several maps from A to B, but by modding out by maps from A 
to A we are left with one map from A to B.) 

For the very special case of a real bipartite t, /i = and attractive interaction there are states 
which do not even break the A to B symmetry. Indeed, according to (3d. 19), there are states for 
which d{x) = and d{x) is constant on the whole lattice. These states are completely translation 
invariant. 

The gauge symmetries of the Hubbard Hamiltonian that we shall consider are the following: 

(ii) Spin SU{2): The action of SU{2) on the Fock space was defined in (3d.40). It is the 
Bogoliubov transform corresponding to a global spin rotation on the one-particle space H. Every 
Hubbard Hamiltonian is invariant under this transformation. 

(iii) Phase U{1): The action of U{1) is given by the Bogoliubov transformation (3d.41). 
Again, every Hubbard Hamiltonian is invariant under this U{1) symmetry. A generator of the U{1) 
transformation (3d. 41) on Fock space is, in fact, the number operator A/". Thus, U{1) invariance of a 
state implies particle conservation, but not necessarily a definite particle number. (See the remark 
after Theorem 2.3.) A normal state is precisely a state where the ^7(1) symmetry is unbroken. 

(iv) Pseudo-spin SU{2): If we are in the real, bipartite case and )U = we have a very large 
symmetry group. In fact the spin SU{2) symmetry is supplemented by the pseudo-spin SU{2) 
symmetry defined in (3d. 50). 

Although the pseudo-spin SU{2) commutes with the spin SU{2), the pseudo-spin is not com- 
pletely independent of the spin. More precisely, the full symmetry group generated by the spin 
SU{2) and pseudo-spin SU{2) transformations is not isomorphic to the group SU{2) x SU{2). 
The full symmetry group is SO {A), see [YZ]. This is so because for the particular matrix w = —1 
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in SU(2), the spin transformation (i.e., the Bogohubov transformation W on Fock space) corre- 
sponding to (w, 1) € SU{2) X SU{2) is identical to the pseudo-spin transformation corresponding 
to {l,w) G SU{2) X SU{2). This easily follows from (3d. 51) with w = —1 since the Bogoliubov 
transformation Ws{w) commutes with the transformation Wbp. This also corresponds with the 
observation that the representations we obtain for the two spins are either both integer or both 
half-integer. For our purposes it is more useful, however, to treat the spin and pseudo-spin trans- 
formations independently. Indeed, there will be cases where the pseudo-spin SU (2) is a broken 
symmetry but the spin SU{2) is not and vice versa. Merely saying that the 5*0(4) symmetry is 
broken would convey much less information. Likewise the U{1) phase symmetry is just a subgroup 
of the pseudo-spin symmetry; the Bogoliubov transforamtion We in (3d.41) is equal to the pseudo- 

spin transformation Wps corresponding to w = ( ^ j in SU{2). We therefore emphasize 

in Table 2 in the case of broken pseudo-spin (t real and bipartite, positive interaction and /x = 0), 
that not only is the pseudo-spin broken but there even exist HF ground states or Gibbs states for 
which the subgroup U{1) is broken. While the pseudo-spin is always broken in this case, there are 
states for which U{1) is unbroken, i.e., there exist normal ground states and Gibbs states. 

(v) Particle- hole Z2: The large symmetry group consisting of the spin SU{2) and pseudo- 
spin SU{2) required reality, bipartiteness and half-filling. If, however, we give up the condition 
of reality we saw in Section Ill.d that the particle-hole symmetry survives as the antiunitary 
transformation given in (3d. 45). 

In case that we have a real t we shall of course not distinguish the unitary and antiunitary 
particle hole transformations as different symmetries. In this case, breaking of pseudo-spin or of 
spin symmetry may or may not imply breaking of particle-hole symmetry (See table 2). According 
to (3d. 47) (with an extra complex conjugation, because we are in the real t case) the Z2 symmetry 
is unbroken if and only if 



i.e., d{x) = and 6{x) is purely imaginary. 

Notice that in the case where the pseudo-spin is a broken symmetry (t real bipartite, attractive 
interaction and fi = 0) the U{1) symmetry must be broken when Z2 is unbroken and conversely 
Z2 must be broken when U{1) is unbroken. Thus, there are two normal states in this case, related 
by a particle-hole transformation, and there is one state (which cannot be a normal state) that is 
invariant under the particle-hole transformation. 

Note that the spatial A-B symmetry (if exists) is broken if and only if the Z2 particle hole 
symmetry is broken. Although, these are different symmetries, in the sense that they act as 
different transformations on the Fock space, they are identical when restricted to the HF ground 
states and Gibbs states. 

The following diagram illustrates the relationships among the four gauge symetries: The 
particle-hole Z2 is a combination of spin and pseudo-spin, while the phase U{1) is really a subgroup 
of the pseudo-spin SU{2). 
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Spin SU{2) Pseudo-spin SU{2) 

\ I \ 

Particle-hole Z2 Phase U{\) 



TABLE NO. 1 Non-Real t 




Symmetries of H 


Broken Symmetries 


Spatial invariance 


Non-bipartite 
and/or 


Spin SU{2), U{1) 


No results except for U = 00, where 
spin SU{2) is broken but U{1) and 
spatial invariance are not. 


Bipartite 
At = 

Attractive interaction* 


Spin SU{2), U{1), 
Z2 


Z2 


A-B symmetry 
broken 


Bipartite 

Ai = 

Repulsive interaction* 


Spin SU{2), U{1), 
Z2 


Spin SU{2), 
Z2 


A-B symmetry 
broken 



* Notice that even though we are in the bipartite case and = 0, the attractive and repulsive 
interaction Hamiltonians are not unitarily equivalent when t is not real 



TABLE NO. 2 Real Bipartite / 




Symmetries of H 


Broken Symmetries 


Spatial invariance 


/' ^ 

At (raclixx- iiiioracljoii 


Spin SU{2). 

U{1) 


U{1) 


Not broken 


/x = 

Attractive interaction* 


Spin SU{2), 
Pseudo spin SU{2) 


Pseudo-spin SU{2) 
(in particular Z2 and 
{/(I) can be broken) 


A-B symmetry 
can be broken 


Repulsive interaction 


Spin SU{2), 
U{1) 


No results except for U = 00, where 
spin SU{2) is broken but U{1) and 
spatial invariance are not. 


/X = 
Repulsive interaction* 


Spin SU{2), 
Pseudo-spin SU{2) 


Spin SU{2) 

(in particular Z2 is 

broken) 


A-B symmetry 
broken 



For real, bipartite t and fj, = the Hamiltonians with attractive and repulsive interactions are 
unitarily equivalent. Because of the (—1)^, a spatial symmetry between the A and B sublattices 
will not commute with this unitary transformation; hence Table 2 fails to be the same for the 
repulsive and attractive cases. 
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TABLE NO. 3 Real, Non-Bipartite t 




Syiii riiclrics of H 


Broken Syiiiriiclrics 


SfHilial irirnruince 


Attractive interaction 
Any /X 


Spin SU{2), 

m 


m 


Not broken 


Repulsive interaction 
Any /i 


Spin SU{2), 

m 


No results except for U = 00, where 
spin SU{2) is broken but U{1) and 
spatial invariance are not. 
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